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Abstract: The gravitino is a promising candidate for cold dark matter. We study cos- 
mological constraints on scenarios in which the gravitino is the lightest supersymmetric 
particle and a charged slepton the next-to-lightest supersymmetric particle (NLSP). We 
obtain new results for the hadronic nucleosynthesis bounds by computing the 4-body de- 
cay of the NLSP slepton into the gravitino, the associated lepton, and a quark-antiquark 
pair. The bounds from the observed dark matter density are refined by taking into account 
gravitinos from both late NLSP decays and thermal scattering in the early Universe. We 
examine the present free-streaming velocity of gravitino dark matter and the limits from 
observations and simulations of cosmic structures. Assuming that the NLSP sleptons freeze 
out with a thermal abundance before their decay, we derive new bounds on the slepton 
and gravitino masses. The implications of the constraints for cosmology and collider phe- 
nomenology are discussed and the potential insights from future experiments are outlined. 
We propose a set of benchmark scenarios with gravitino dark matter and long-lived charged 
NLSP sleptons and describe prospects for the Large Hadron Collider and the International 
Linear Collider. 
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1. Introduction 

Supersymmetric extensions of the standard model of particle physics have remarkable prop- 
erties [1, 2]. In particular, a theory invariant under general coordinate transformations, 
i.e., supergravity (SUGRA), is obtained once supersymmetry (SUSY) is promoted from a 
global to a local symmetry. 
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The gravitino is the gauge field of local SUSY transformations. It is the spin-3/2 
superpartner of the graviton and, thus, a singlet with respect to the gauge groups of 
the standard model. For unbroken SUSY, the gravitino is massless and interacts only 
gravitationally. Once local SUSY is broken spontaneously, the gravitino acquires a mass 
and the interactions of the spin- 1/2 goldstino (which is the Goldstone spinor associated 
with the breaking of global SUSY) through the super-Higgs mechanism. The gravitino 
mass depends strongly on the SUSY-breaking scheme and can range from the eV scale to 
scales beyond the TeV region [1, 3, 4, 5]. For example, the gravitino mass is typically less 
than 100 MeV in gauge-mediated SUSY breaking schemes [3] and is expected to be in the 
GeV to TeV range in gravity-mediated schemes [1]. The gravitino interactions — given by 
the SUGRA Lagrangian [6, 2] — are suppressed by inverse powers of the (reduced) Planck 
scale [7]. At energies well above the gravitino mass, these interactions are enhanced as the 
spin-1/2 goldstino components, which enter through the super-Higgs mechanism, become 
dominant. 

We assume in the following that the gravitino is the lightest supersymmetric particle 
(LSP) which is stable due to i?-parity conservation. If heavier than about 100 keV, the 
gravitino LSP is a natural candidate for cold dark matter (see [8, 9, 10, 11, 12, 13, 14, 15, 
16, 17, 18, 19, 20, 21, 22, 23, 24, 25] and references therein). Gravitino dark matter can 
originate from both thermal production in the very early Universe [11, 12, 14, 16] and late 
decays of the next-to-lightest supersymmetric particle (NLSP) [13, 15, 17, 18, 21, 22]. We 
do not consider other more model-dependent gravitino sources such as inflaton decay [26]. 

The gravitino LSP scenario is subject to severe cosmological constraints. Upper limits 
on the reheating temperature after inflation can be obtained since the relic abundance 
of thermally produced gravitinos cannot exceed the observed abundance of dark mat- 
ter [11, 15, 16, 23, 25]. Likewise, one can derive upper limits on the gravitino and NLSP 
masses since the gravitino density from NLSP decays is also bounded from above by the ob- 
served dark matter density [15, 22]. Moreover, there are limits on the present free-streaming 
velocity of the gravitino dark matter from observations and simulations of cosmic struc- 
tures. Depending on the origin of the gravitinos, these limits can further restrict the allowed 
range of the gravitino mass and the mass of the NLSP [13, 27]. 

In addition to the limits from the properties of dark matter, there are also constraints 
from the observed abundances of the primordial light elements (D, He, Li) [20, 21, 22, 23]. 
Because of the extremely weak interactions of the gravitino, the NLSP typically has a long 
lifetime before it decays into the gravitino. If these decays occur during or after primordial 
nucleosynthesis, the standard model particles emitted in addition to the gravitino can affect 
the abundance of the primordial light elements. The effects of both electromagnetic and 
hadronic showers can be important depending on the lifetime of the NLSP (see [28, 29] 
and references therein). Bounds on the gravitino and NLSP masses can be derived by 
demanding that the successful predictions of primordial nucleosynthesis be preserved. 

Nucleosynthesis constraints have been examined for neutralino NLSP and slepton 
NLSP scenarios. While most of the existing investigations focus on the effects of elec- 
tromagnetic energy injection (see [20, 30, 31] and references therein), the bounds from 
hadronic energy injection can be much more severe [28, 29, 32, 33, 34, 35]. In particular, 
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the hadronic constraints strongly restrict the mass spectrum in the bino-like neutralino 
NLSP case [22, 23]. For the slepton NLSP case, the hadronic constraints have also been 
estimated and found to be much weaker but still significant in much of the parameter 
space [21, 22, 23]. 

In this paper we study systematically the cosmological constraints on the gravitino LSP 
scenario for the case of a charged slepton NLSP. Indeed, the lighter stau t\ can appear nat- 
urally as the lightest standard model superpartner in frameworks in which the spectrum at 
low energies is derived from SUGRA predictions at the scale of grand unification. Scenarios 
with a gravitino LSP and a long-lived charged slepton NLSP are also particularly promising 
for phenomenology at future colliders [36, 37, 38, 39, 40, 41, 42, 22, 43, 44, 45, 25, 46]. 

We present new results for the hadronic energy release in late decays of a charged 
slepton NLSP. We obtain these results from an exact computation of the 4-body decay 
of the slepton NLSP into the gravitino, the corresponding lepton, and a quark-antiquark 
pair. The corrections with respect to the previous estimate [21] are indicated explicitly. 

Based on our new results, we derive upper limits on the abundance of the charged NLSP 
slepton before its decay. These limits lead to new bounds on the NLSP and gravitino 
masses in scenarios in which the NLSP sleptons freeze out with a thermal abundance 
before their decay. We also refine the bounds from the observed abundance of dark matter 
by taking into account gravitinos from both slepton NLSP decays [13, 15, 17, 18, 21, 
22] and thermal scattering in the early Universe [11, 12, 14, 16]. Upper limits on the 
reheating temperature after inflation are extracted based on the gauge-invariant treatment 
of thermal gravitino production [16]. For gravitinos from both NLSP decays and thermal 
production, we compute the present free-streaming velocity and confront it with limits from 
observations and simulations of cosmic structures. For each of the considered constraints, 
we explore associated uncertainties and outline potential refinements from astrophysical 
observations and experiments at future colliders. 

We discuss implications of the cosmological constraints for gravitino dark matter, col- 
lider phenomenology, the SUSY breaking scale, and the generation of the baryon asymme- 
try in thermal leptogenesis [47, 48]. Based on our findings, we propose and examine a set 
of benchmark scenarios with gravitino dark matter and long-lived charged NLSP sleptons. 
We consider a wide range of gravitino, slepton NLSP, and gluino masses, and our study is 
not restricted to a constrained framework such as the constrained minimal super symmetric 
extension of the standard model (CMSSM) [20, 23]. 

The rest of this paper is organized as follows. In Sec. 2 we consider the 2-body decay 
of the charged slepton NLSP into the gravitino and the associated lepton. This decay 
mode governs the lifetime of the NLSP and the bounds from late electromagnetic energy 
injection. In Sec. 3 we present our results on the hadronic 4-body decay of the slepton 
NLSP. Section 4 reviews the nucleosynthesis constraints for late decaying particles and 
provides our new upper limits on the abundance of the charged slepton NLSP before 
its decay. In Sec. 5 we study the relic density and the present free-streaming velocity of 
gravitinos from NLSP decays and the corresponding bounds. An update of the cosmological 
constraints on the gravitino and slepton NLSP masses is given in Sec. 6. In Sec. 7 we discuss 
the thermal production of gravitinos, their present free-streaming velocity, and the bounds 



-3- 



on the reheating temperature after inflation. Section 8 describes the possible interplay 
between astrophysics and collider phenomenology in view of the cosmological constraints. 
We summarize the implications for cosmology and collider phenomenology by proposing ten 
viable benchmark scenarios with gravitino dark matter and a long-lived charged slepton. 

2. 2-Body Decays of Charged NLSP Sleptons 

In this section we consider the 2-body decay of a charged NLSP slepton into the corre- 
sponding lepton and the gravitino LSP. We study the lifetime of the slepton NLSP, which 
is governed by this decay mode, and discuss the release of electromagnetic and hadronic 
energy resulting from such 2-body decays in view of the nucleosynthesis constraints. 

We concentrate on the scenario where the lighter stau t\ is the NLSP, but address 
also the alternative selectron or smuon NLSP cases. Indeed, the mass spectrum of the 
superparticles could be such that the masses of the lighter selectron ei, the lighter smuon 
Jli, and the lighter stau t\ are nearly degenerate, nif « m^p,- Then, the decays of these 
selectrons, smuons, and staus can be equally important for the constraints from primordial 
nucleosynthesis. However, if the stau is significantly lighter than its first- and second- 
generation counterparts, the effects of the selectron and smuon decays on the primordial 
light elements will be negligible. This results from negligible abundances of selectrons and 
smuons at the time when the stau NLSP decays into the gravitino LSP; cf. Sec. 6.1. 

The two staus are in general a linear combination of tr and tl, which are the superpart- 
ners of the right-handed and left-handed tau lepton, respectively: = cos # r ?R-|-sin# r rL. 
We focus for simplicity on a purely 'right-handed' lighter stau, t\ = tr. This is a good 
approximation at least for small values of tan/3 = v u /vd, where v u d = denote the 

vacuum expectation values of the two neutral Higgs scalar fields of the minimal supersym- 
metric standard model (MSSM). Then, the neutralino-stau coupling is dominated by the 
bino coupling. We also assume for simplicity that the lightest neutralino is a pure bino. 

2.1 Lifetime of the Charged Slepton NLSP 

In the gravitino LSP scenario, the 2-body decay r — > tG is the main decay mode of the 
stau NLSP. Neglecting the r mass, the SUGRA Lagrangian [6, 2] leads to the following 
tree-level result for the decay rate 



which is also valid for any mixture r = cos 6 t tr + sin 9 t ti,. To get expression (2.2), we have 
used the value of the reduced Planck mass M P1 = (8^G N r 1/2 = 2.435 x 10 18 GeV as ob- 
tained from macroscopic measurements of Newton's constant [7] Gn = 6.709x 10~ 39 GeV -2 . 




(2.1) 



(2.2) 



Thus, the stau NLSP lifetime, Tf « 1/T(tr — > tG), is determined by the gravitino mass 
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Stau NLSP Lifetime 




Figure 1: The lifetime of the stau NLSP as a function of its mass m,f for values of the gravitino 
mass ttiq ranging from 10 MeV up to 1 TeV. 

and the stau NLSP mass mf. With the replacement mf — > m^/x, the results for the 
selectron NLSP and smuon NLSP cases are identical to the one given above. 

In Fig. 1 we plot the lifetime of the stau NLSP as a function of its mass mf for values 
of the gravitino mass ranging from 10 MeV up to 1 TeV. This figure shows the strong 
sensitivity of the stau NLSP lifetime on and mf. Indeed, Tf can be in the range 
from 0.01s to 10 12 s and beyond. The contribution from the longitudinal helicity-1/2 (or 
goldstino) components to the decay rate dominates for mf > 3mg. For values of mf close 
to mg, these contributions become less important and the phase space shrinks so that the 
stau NLSP lifetime tends to infinity as illustrated for mg = 100 GeV and mg = 1 TeV. 

2.2 Electromagnetic and Hadronic Energy Release from the 2-Body NLSP De- 
cay 

Depending on the lifetime Tf, the tau emitted in the stau NLSP decay can affect the 
abundance of the light elements produced in primordial nucleosynthesis. At the time of its 
decay, the stau NLSP is non-relativistic and can be considered to be at rest with respect 
to the background. In the NLSP rest frame or comoving frame, the initial energy of the 
tau (with a mass of m T = 1.777 GeV) from the 2-body decay is given by 

m~ — m?~ + ml 



Accordingly, the emitted tau has a time-dilated lifetime of r r = 2.9 x 10 13 (E T /m T ) s . 
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If the tau scatters off the background plasma before its decay, the 2-body decay of one 
stau NLSP gives no hadronic energy release, ehadC? — ► tG) = 0, but an electromagnetic 
energy release of e cm (r tG) = E T , i.e., the full initial tau energy (2.3) ends up in an 
electromagnetic shower. However, when the temperature of the Universe drops below about 
0.5 MeV, which corresponds to a cosmic time of a few seconds, the interaction time of the 
tau with the background plasma exceeds its time-dilated lifetime even for E T = 0(1 TeV). 
One then has to consider the tau decays and the effects of the associated decay products. 

In each tau decay, at least one neutrino is emitted. Since neutrinos are only weakly 
interacting, the effect of their injection is typically subleading [49, 50]. With (basically 
invisible) energy carried away by neutrinos, the electromagnetic energy release e cm (r — ► tG) 
can be reduced down to about one third of the initial tau energy (2.3) [18]. 

The tau has a sizeable number of decay modes. In addition to neutrinos, one finds 
that also electrons and unstable particles such as fi, tt^, 7r°, K ± , and K° appear as decay 
products. For the computation of e emi had('r -»tG), one therefore has to examine not only 
the interactions of these secondary particles but also the decays of the unstable particles 
and the interactions of their decay products with the primordial plasma. Depending on 
the envisaged level of rigor, this can be a tedious task. 

Let us here describe the effects of the mesons emitted in the main decay modes of 
the tau lepton without going into the details. The neutral 7r° meson contributes only to 
electromagnetic energy injection because it decays purely electromagnetically before inter- 
acting with the background plasma. The charged mesons ir^ and also contribute to 
electromagnetic energy injection. Indeed, their electromagnetic interaction time is signifi- 
cantly shorter than their hadronic one [32]. These mesons thus deposit basically all of their 
kinetic energy in the electromagnetic cascades when scattering off the background plasma. 
After these mesons are stopped, they still can induce proton-neutron interconversion pro- 
cesses, which can affect the abundances of the primordial light elements. Also, the neutral 
K® mesons, which cannot be stopped electromagnetically, can trigger proton-neutron in- 
terconversion processes. At cosmic times of t > 100 s or, equivalently, temperatures of 
T < 0.1 MeV, however, the emitted 7r ± , K ± , and K° mesons typically decay before inter- 
acting hadronically [29] . Accordingly, the hadronic effects of the mesons from the tau decays 
can only be relevant for stau NLSP decays at temperatures of 0.1 MeV < T < 0.5 MeV, 
which corresponds to cosmic times of 3 s < t < 100 s. 

In this paper we will concentrate on decays at late times of t > 100 s, where the 
hadronic effects of the emitted mesons can be neglected so that ehadC? -> tG) = 0. Then, 
the 4-body decay of the stau NLSP into the gravitino, the tau, and a quark-antiquark pair 
governs the constraints from late hadronic energy injection. The hadronic energy release 
from this 4-body decay will be computed in the next section. 

Upper limits on electromagnetic energy release become most severe at cosmic times 
of 10 7 s < t < 10 12 s, as will be shown explicitly in Sec. 4.2. The resulting constraints 
are governed by the 2-body decay r — > tG and were studied for 6 e m 

(r — > r G) = 0.5 .Et- 
ui [18, 21, 22, 23] and for e cm (r -trG) = 0.3 E T in [20]. The true value of e cm (r -4rG) 
can deviate from these values depending on the lifetime of the stau NLSP. We will explore 
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Figure 2: The electromagnetic energy release from the 2-body decay of a single stau NLSP as a 
function of m,f for mg < 1 GeV and mg = 100 GcV, 300 GeV, and 1 TeV. The solid, dashed, and 
dotted lines are obtained with e em (T — > tG) = E T , 0.5 E T , and 0.3 E T , respectively 



the electromagnetic constraints again only for representative values: 

e cm (r ^ tG) =xE T with x = 0.3, 0.5, 1 . (2.4) 

In Fig. 2 we show e cm (r — > r G) = E T (solid lines), 0.5 E T (dashed lines), and 0.3 E T (dotted 
lines) as a function of m,f for < 1 GeV and mg = 100 GeV, 300 GeV, and 1 TeV. 
Note that the true electromagnetic energy release is always below e cm (r — ► tG) = E T for a 
stau NLSP that decays at a cosmic time of 100 s < t < 10 12 s. Accordingly, the constraints 
obtained for e em (r — ► r G) = E T are overly restrictive for stau NLSP scenarios. 

In the selectron NLSP case, the electromagnetic energy injection is governed by the 
2-body decay e — > e G in which a (stable) electron is emitted with an energy of E e = (m| — 
m% + m 2 )/(2 771,5). Accordingly, an electromagnetic and hadronic energy of respectively 

£em (e — ► e G) ~ E e and ehad (e — > e G) ~ is released in the 2-body decay of one selectron 
NLSP. For the case of the selectron NLSP, the electromagnetic energy release is described 
by the same solid curves in Fig. 2 since the differences resulting from the different lepton 
masses are negligible. 

In the smuon NLSP case, the 2-body decay Jl — ► fxG injects a muon with an energy of 
E^ = (m~ — m 2 ~ + m^)/(2mp) into the primordial plasma. This muon has a time-dilated 
lifetime of 2.2 x 10" 6 {E^/m^) s . As long as the temperature of the Universe is above 
about 0.3 keV, which corresponds to cosmic times t < 1.5 x 10 4 s, this muon initiates an 
electromagnetic cascade so that e era (jl — ► fj,G) ~ while ehad(A* — ► // G) — . If the muon 
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Figure 3: The 4-body NLSP decay tr -> rGqq 

is emitted at later times or, equivalently, lower temperatures, it decays via (i — ► ev e v^ before 
interacting with the background plasma. Accordingly, the electromagnetic energy release 
is reduced to about e cm (Jl — ► jj,G) f» 0.3 £^ for smuon NLSP decays at i > 1.5 x 10 4 s. 
The hadronic energy release remains at ehad(^ — ► /■*(?) ~ 0. For such late decays, the 
electromagnetic energy release as a function of the smuon mass can be read from the 
dotted curves in Fig. 2 since the effect of the different lepton masses is again negligible. 

3. 4-Body Decays of Charged NLSP Sleptons into Hadrons 

In this section we compute the 4-body decay of a charged NLSP slepton into the corre- 
sponding lepton, the gravitino LSP, and a quark-antiquark pair. We present our results for 
the branching ratio of this decay mode, the associated hadronic energy spectrum, and the 
resulting hadronic energy release per slepton NLSP decay. We compare the results from 
our exact treatment of the 4-body decay with the previous estimate derived in Ref. [21]. 

We compute the leading contribution to the 4-body decay 1r — > lGqq. To be specific, 
we discuss our results for the scenario in which the tr is the NLSP. Nevertheless, the results 
given in this section can be applied directly to eR or ptR NLSP scenarios by performing 
obvious substitutions such as nif — > mg,^. 

The Feynman diagrams of the leading contributions to tr — > rGqq are shown in Fig. 3. 
The couplings of the gravitino to the stau, tau, bino (£>), and the electroweak gauge bosons 
are given by the SUGRA Lagrangian [6, 2]. In our computation of the squared matrix 
element, the exchange of virtual photons (7*) and virtual Z bosons (Z*) is included. The 
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width of the Z boson is taken into account by using the Breit-Wigner form of the Z-boson 
propagator. We neglect the lepton mass. We also assume that the lightest neutralino is a 
pure bino as already mentioned. The other neutralinos and the squarks are assumed to be 
much heavier than the stau NLSP and the bino. Depending on the invariant mass of the 
quark-antiquark pair m qq , up to five quark flavors are considered. The contributions from 
tt-final states, which can appear for m qq > 350 GeV, are subleading. The masses of the 
quarks are taken into account only in the phase space integration. Only quark-antiquark 
pairs with an invariant mass of at least = 2 GeV, i.e. the mass of a nucleon pair, are 
considered. This choice for the cut on m qq will be motivated below. Performing the phase 
space integrations numerically, we obtain the results presented in Figs. 4 to 8 below. 
The branching ratio of the 4-body decay tr — > rGqq with rn qq > is given by 



Tfrp ->rGqq; m c r{r . 
BR(tr -rGqq;m;f)E qq ' , (3.1) 

TR 

where the total width of the stau NLSP r~° tal is dominated by the 2-body decay 

ptotal ~ r (f B T G) . (3.2) 
The partial width of the 4-body decay with a cut on the invariant mass m qq is defined as 

r(r R -> rGqq;m^ 1 = / dm m . (3.3) 



'qq 



In Fig. 4 the branching ratio of the decay tr — > rGqq is shown as a function of m,f for 
m G ~ ^ GeV (left) and = 100 GeV (right). The curves are obtained with m^* = 2 GeV 
and a bino mass of = 1.1 mf. The contributions from pure 7* and pure Z* exchange 
are indicated by the dashed and dotted lines, respectively. The solid lines present the full 
results in which also the j*-Z* interference is taken into account. For m.f — > mg, the 
phase space shrinks so that the considered branching ratio tends to zero as illustrated for 
mg = 100 GeV. 

In Fig. 5 we show the branching ratio of tr — > rGqq as a function of mf for values of 
up to 1 TeV, where, again, rrf^ = 2 GeV and = 1.1 mf. The full results are given 

by the solid lines and the respective contributions from pure Z* exchange by the dotted 

lines. 

In previous studies, the hadronic branching ratio was estimated from a real Z boson 
in the final state and its hadronic branching fraction [21]. These estimates (obtained in 
the zero- width approximation) differ very little from the dotted lines shown in Figs. 4 
and 5. In other words, the hadronic branching ratio was significantly underestimated 
in Refs. [21, 22, 23] towards smaller values of mf , where the photon exchange gives the 
dominant contribution. 

The energy spectrum of the emitted hadrons is also important for the investigation of 
the nucleosynthesis constraints. As mentioned in the previous section, we concentrate on 
decays at cosmic times of t > 100 s, where the hadronic effects of the mesons from the tau 
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m f [TeV] ra f [TeV] 

Figure 4: The branching ratio of the 4-body decay tr — > rGqq as a function of rrif for < 1 GeV 
(left) and mg = 100 GeV (right). The solid lines give the full results including j*-Z* interfer- 
ence. The dashed and dotted lines show the contributions from pure 7* and pure Z* exchange, 
respectively. The curves are obtained with m^* = 2 GeV and a bino mass of = 1.1 m?. For 
comparison, the previous estimates of the branching ratio inferred from a real Z boson (i.e., in the 
zero width approximation) alone [21] coincide basically with the dotted lines. 



decays can be neglected. Thus, the initial hadronic energy is determined by the invariant 
mass m q q of the quark-antiquark pair emitted in the considered 4-body decay. 

In Fig. 6 the energy spectrum of the quark-antiquark pair normalized to the partial 
width of the 4-body decay, 

1 dTfa^rGqq) ^ 



r(7k->TGqq;m^ = 2 GeV) dm m 

is shown for (m d ,m f ) = (1 GeV, 1 TeV), (1 GeV, 0.2 TeV), (100 GeV, 0.2 TeV), and 
(10 GeV, 0.1 TeV) by the solid, dashed, dotted, and dash-dotted curves, respectively, where 
wig = 1.1 mf. The solid, dashed, and dotted curves illustrate that the importance of the 
Z-boson resonance at m q q = Mi grows with increasing mass difference Am = mf — mg. 
For a mass difference Am less than the Z-boson mass, however, the Z boson cannot become 
resonant and the photon exchange dominates as illustrated by the dash-dotted curve. The 
kinks at m qq ~ 4 GeV and 10 GeV mark the threshold energies for the production of cc 
and bb pairs respectively. 

Having computed the initial energy spectrum of the quark-antiquark pairs for given 
values of mf and mg, one should consider the fragmentation of the initial quark-antiquark 
pairs into hadrons. The resulting meson and nucleon spectra should then be taken into 
account when computing the abundances of the primordial light elements; cf. [29, 51]. We 



- 10 - 



Branching Ratio of tr — > rGqq 

0.01 p 1 — |— 



t 



10" 



10" 



10 




rn n [TeV] 



m~ — 1.1 m,T 



i i I I i i i 



0.1 0.2 0.3 0.5 1 2 3 
m f [TeV] 



10 



Figure 5: The branching ratio of the 4-body decay tr — > rGqq as a function of mf for < 1 GeV 
and uIq — 0.1, 0.2, 0.3, 0.5, 0.7, and 1 TeV (from the left to the right). The solid lines present the 
full results and the dotted lines the contributions from pure Z* exchange for comparison with the 
previous estimates. The curves are obtained with m™' = 2 GeV and a bino mass of nig = 1.1 mf. 



will not follow this procedure. Instead, we will adopt the limits on the initial hadronic 
energy release from [28, 29], where the hadronization of primary partons is already taken 
into account. These limits are most severe for decays at cosmic times t > 100 s. For such 
late decays, the emitted nucleons govern the constaints since — as already mentioned — 
the mesons typically decay before interacting with the background nuclei. We therefore 
consider only quark-antiquark pairs with an invariant mass greater than the mass of a pair 
of nucleons, m%j£ = 2 GeV. 

Our study of the hadronic nucleosynthesis constraints is based on the computation of 
the average initial hadronic energy release in the 4-body decay of a single stau NLSP 

Chad (tr r G q q) = dm m m q q ^— . (3.5) 

In Fig. 7 this quantity is shown as a function of mf for < 1 GeV (left) and mg = 
100 GeV (right), where m^' = 2 GeV and = 1.1 mf. The solid lines present the 
full results and the dashed and dotted lines the contributions from pure 7* and pure Z* 
exchange, respectively. In Fig. 8, the solid lines show the full results for gravitino masses 
up to 1 TeV. In addition, the dotted lines show the curves obtained with the estimate from 
Ref. [21] 

e h ad(FST '04) = - (m f - m^) BR(tr rGZ)BR(Z qq) , (3.6) 
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Figure 6: The normalized energy spectrum of quark-antiquark pairs emitted in the 4-body decay 
tr — ► rGqq with an invariant mass m q q. The solid, dashed, dotted, and dash-dotted curves are 
obtained respectively for (mg,m f ) = (1 GeV, 1 TeV), (1 GeV, 0.2 TeV), (100 GeV, 0.2 TeV), and 
(10 GeV, 0.1 TeV), where m s = l.lm f . 



where the energy of the quark-antiquark pair was approximated as [m,f — m£,)/3 and 
only Z-boson exchange in the zero- width approximation was taken into account. From the 
comparison, we find that the hadronic energy release has been underestimated towards 
smaller values of nif due to the missing contribution from photon exchange and overes- 
timated for nif 3> mg due to the monochromatic approximation of the energy of the 
quark-antiquark pair. Accordingly, one has to update the hadronic constraints inferred 
from the estimate (3.6) such as the ones given in [21, 22, 23]. 



4. Constraints from Primordial Nucleosynthesis 

In this section we review the nucleosynthesis constraints on hadronic and electromagnetic 
energy release from late NLSP decays. We use these constraints to derive upper limits on 
the abundance of charged slepton NLSPs before their decay. 

4.1 Primordial Nucleosynthesis and Effects of Late NLSP Decays 

Primordial nucleosynthesis, or big-bang nucleosynthesis (BBN), is one of the cornerstones 
of modern cosmology. With the observed abundance of D and 4 He in the Universe, the 
standard BBN scenario provides an estimation of the baryon-to-photon ratio, rj = ub/tIj, 
which agrees with that inferred from studies of the cosmic microwave background (CMB). 
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Figure 7: The average hadronic energy release from the 4-body decay of a single stau NLSP as 
a function of mf for mg < 1 GeV (left) and mg = 100 GeV (right), where m^* = 2 GeV and 
vrig — 1.1 mf. The solid lines show the full results, and the dashed and dotted lines the contributions 
from pure 7* and pure Z* exchange, respectively. 

As mentioned already in the Introduction, this concordance imposes constraints on the 
considered gravitino LSP scenario since the electromagnetic and hadronic energy released 
in late NLSP decays can modify the standard BBN scenario substantially. 

The dominant mechanism affecting the standard BBN scenario depends on the time 
t at which the electromagnetic or hadronic energy is injected. In this study we use the 
sudden decay approximation by assuming that the release of electromagnetic and hadronic 
energy from the NLSP decay occurs when the cosmic time t equals the NLSP lifetime 
tnlsp- For Is < tnlsp ^5 100s, energetic hadrons are stopped efficiently through elec- 
tromagnetic interactions so that the direct destruction of light elements is subdominant. 
The presence of additional slow hadrons can still change the ratio of protons to neutrons 
through interconversion processes and thus affect the abundance of the light elements. For 
100 s < tnlsp ^5 10 7 s, energetic hadrons and, in particular, neutrons cannot be slowed 
down significantly. Accordingly, they can reprocess efficiently the produced light elements 
through hadro-dissociation processes. The effect of electromagnetic energy release is negli- 
gible for tnlsp ^ 10 4 s as the interaction with the background particles thermalizes quickly 
any high-energy photons or leptons emitted in the NLSP decay. Towards later times, elec- 
tromagnetic energy release becomes important. For 10 7 s < tnlsp ^5 10 12 s, the reprocessing 
of light elements through energetic electromagnetic showers, i.e., photo-dissociation, can 
become more significant than hadro-dissociation. (For details, see [32, 29] and references 
therein.) 
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Figure 8: The average hadronic energy release from the 4-body decay of a single stau NLSP 
as a function of mj for = 0.1, 0.2, 0.3, 0.5, 0.7, and 1 TeV (solid lines from the left to the 
right), where = 2 GeV and trig = 1.1 nif . The corresponding previous estimates from [21] are 
indicated by the dotted lines (FST '04). 



The constraints on electromagnetic and hadronic energy release can be quantified con- 
veniently in terms of upper bounds on the quantities 

£em,had = €em,had Ynlsp • (4.1) 

Here e em ,had is the (average) electromagnetic/hadronic energy emitted in a single NLSP 
decay studied for charged slepton NLSP scenarios in the previous sections. The NLSP 
yield prior to decay Ynlsp is obtained by dividing the NLSP number density unlsp prior 
to decay with the total entropy density of the Universe s: 

Ynlsp = ^ • (4.2) 
s 

In this way, one scales out the effect of the expansion of the Universe on the NLSP abun- 
dance. 

We concentrate on the bounds from the primordial abundances of D and 4 He as they 
are the most reliable ones. The bounds from 3 He, 6 Li, and 7 Li could be even more severe. 
However, the observed abundances of these elements are still subject to serious systematic 
uncertainties. 

We work with the recent numerical constraints on electromagnetic and hadronic en- 
ergy release given in Refs. [28, 29]. These constraints were extracted with the following 
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observational values of the light element abundances (Id error): 



(n D /n H )mean = (2.78±g;^) x 10" 5 , 
(n D /n H )high = (3.98±°;|?) x 10~ 5 , 
y p (FO) = 0.238 ± (0.002) stat ± (0.005) syst , 
F P (IT) = 0.242 ± (0.002) stat [±(0.005) syst ] 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



Here «x denotes the primordial number density of element X and Y p the primordial value 
of the mass fraction of 4 He. Observe that two values for wd/^h were considered. The one 
given in (4.3) corresponds to the mean of the existing observations and the other given 
in (4.4) to the highest value among these observations. The latter of these was used to 
derive a conservative constraint; cf. Sec. IIA of Ref. [29] and references therein. For Y p , 
the values from the analysis of Fields and Olive (FO) [53] and from the one of Izotov and 
Thuan (IT) [54] were used due to the sizeable difference. 

We also consider the recent numerical constraints on electromagnetic energy release ex- 
tracted from more conservative observational bounds on the light element abundances [31] 



For nj)/nii, the upper limit is the 2a upper limit to the highest value reliably observed. 
The lower limit is adopted from the observed present abundance in the interstellar medium 
as galactic processes can only destroy deuterium. For the 4 He abundance, the range is 
obtained from the analysis of Fields and Olive [53] by taking the 2a range with errors 
added in quadrature. 

In Fig. 9 we show the numerically obtained upper limits on hadronic and electromag- 
netic energy release as a function of the NLSP lifetime tnlsp • The upper part of the figure 
shows exclusion limits on hadronic energy release £had = e had^NLSP taken from Fig. 39 
of Ref. [29]. The solid (black), dashed (red), dotted (blue), and thin dash-dotted (yel- 
low) lines represent the limits inferred from the observational bounds (4.3), (4.5), (4.6), 
and (4.4), respectively. These limits were obtained at the 95% confidence level (CL) with a 
baryon asymmetry of rj = (6.1 ± 0.3) x 10~ 10 by considering a particle of lifetime tnlsp an d 
mass ITeV decaying with a hadronic branching ratio of Shad = 1- Note that these limits 
were found to be rather insensitive to the mass of the decaying particle. The lower part 
of Fig. 9 shows exclusion limits on electromagnetic energy release ^ em = e cm 

Fnlsp- The 

solid (black), dashed (red), and dotted (blue) lines are inferred respectively from the ob- 
servational bounds (4.3), (4.5), and (4.6). These curves are taken from Fig. 42 of Ref. [29]. 
They indicate the limits at the 95% CL obtained with rj = (6.1 ±0.3) x 10~ 10 for a particle 
of lifetime tnlsp and mass 1 TeV decaying purely electromagnetically, B em = 1. In ad- 
dition, the exclusion limits on £em from the more conservative observational bounds (4.7) 
and (4.8) are indicated by the thin dash-dotted (yellow) and dotted (yellow) lines, respec- 
tively Accounting for the different normalizations, we adapt these curves (obtained with 
i| = 6x 10~ 10 ) from Figs. 6 and 7 of Ref. [31]. 



1.3 x 10" 5 < n D /n H < 5.3 x 10" 
0.227 < Y p (< 0.249) . 



-5 



(4.7) 
(4.8) 
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Upper Limits on Hadronic and Electromagnetic Energy Release 
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Figure 9: Upper limits on hadronic and electromagnetic energy release from NLSP decays as a 
function of the NLSP lifetime tnlsp • The upper plot shows the bounds on ehad ^nlsp adapted from 
Fig. 39 of Ref. [29] and the lower plot the bounds on e em Inlsp adapted from Fig. 42 of Ref. [29] 
and Figs. 6 and 7 of Ref. [31]. The solid (black), dashed (red), and dotted (blue) lines result from 
the observational constraints (4.3), (4.5), and (4.6), respectively. The thin dash-dotted (yellow) 
lines in the upper and lower plot are obtained from (4.4) and (4.7), respectively. The thin dotted 
(yellow) curve represents the limit from (4.8). 

The upper part of Fig. 9 shows that the limits on hadronic energy release become most 
severe for tnlsp ^ 100 s, where hadro-dissociation processes govern the constraints. As 
mentioned in the previous sections, we concentrate on this region. We do not consider the 
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Figure 10: Contours of the hadronic energy ehad (solid lines, left plot) and the electromagnetic 
energy e cm = 0.5 E T (solid lines, right plot) released in the decay of a single right-handed stau in the 
parameter space spanned by and mf. The dashed contours indicate stau lifetimes of t? = 0.01, 
1, 100, 10 4 , and 10 7 s (from the upper left to the lower right). 

inter conversion processes triggered by mesons, which become relevant for tnlsp i$ 100s. 
For 10 7 s < tnlsp i$ 10 12 s, the upper limits on purely electromagnetic energy release 
become as severe as the hadronic ones as can be seen from the lower part of Fig. 9. For 
more details on these limits and their interpretation, we refer to the original papers [29, 31]. 

In the following we will work with one severe and one conservative upper limit for 
both hadronic and electromagnetic energy release irrespective of the details. The severe 
limits are given by the solid (black) lines in conjunction with the dashed (red) lines. The 
conservative limits are obtained from combinations of the thin dash-dotted (yellow) and 
dotted (blue) curves. 

4.2 Upper Limits on the Abundance of Charged Slepton NLSPs before their 
Decay 

Let us focus again on scenarios with a charged slepton NLSP decaying late into the gravitino 
LSP. The corresponding (average) electromagnetic/hadronic energy emitted in a single 
charged slepton NLSP decay e em! had was studied in Sees. 2 and 3. From the nucleosynthesis 
bounds on £ e m,had(= Cem,had ^nlsp) discussed above, we now derive upper limits on the 
yield of charged slepton NLSPs before their decay. Once the yield is computed, the limits 
presented in this section can be transformed into bounds on the masses of the charged 
slepton NLSP and the gravitino LSP as will be demonstrated in Sec. 6. Since the extracted 
limits are sensitive to the nucleosynthesis bounds, we will show our results for both the 
severe and the conservative limits on £ e m,had defined above. 

For a better understanding of the limits given below, Fig. 10 shows contours of the life- 
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Upper Limits on Y-f- — / s prior to Stau NLSP Decays 
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Figure 11: Upper limits on the yield of right-handed stau NLSPs before their decay, Y~^* d = 
£had X / e had, as obtained from the severe (left) and the conservative (right) upper limits on late 
hadronic energy injection specified at the end of Sec. 4.1. The contour lines indicate Y-^j, = 
3 x 10~ 13 , 10~ 12 , 10~ n , and 1CT 10 (from the upper right to the lower left). A darker shading 
implies a smaller value of Y~^ d . 



time of the right-handed stau NLSP (dashed lines) together with contours of the hadronic 
energy release ehad (solid lines, left) obtained from (3.5) and the representative electromag- 
netic energy release e em = 0.5i? T (solid lines, right) in the parameter space spanned by mg 
and mf. At each point in the considered parameter space, one finds ehad <C ^cm and might 
naively think that the constraints from hadronic energy injection are negligible. This is 
not the case since Chief ^ f° r T f ~ 10 7 s - (For a gi ven T f = Tnlsp> the upper limits 
on Chad and £em can be read directly from the curves shown in the upper and lower part of 
Fig. 9 respectively.) 

Let us comment on the region in the (mg, mf ) plane for which constraints will be 
provided. Studying gravitino LSP scenarios, we are interested only in the region < mf. 
Moreover, only stau masses mf > 100 GeV are considered since the searches for long-lived 
charged particles at LEP show that masses less than 97.5 GeV can be excluded at 95% CL 
for smuons and staus with a lifetime tj > 1CP 6 s [55]. (The corresponding mass bound for 
long-lived selectrons is 91 GeV.) As already mentioned, we concentrate on late decays, 
Tf > 100 s, for which the hadronic constraints are most severe due to hadro-dissociation 
processes. No limits will be provided for stau lifetimes of Tf < 100 s, where the constraints 
on Chad are weaker anyhow as can be seen in the upper part of Fig. 9. 

In Fig. 11 we give the upper limits on the yield of right-handed stau NLSPs before 
their decay, Y-^^ = Chad X / e had ; as obtained from the severe (left) and the conservative 
(right) upper limit on late hadronic energy injection specified at the end of Sec. 4.1. With 
e had computed from (3.5), this is one of the main results of this paper. We show contours 
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for the values l-f^d = 3 x 1CT 13 , 1CT 12 , 10~ n , and 10~ 10 (from the upper right to the 
lower left), where a darker shading implies a smaller value of Y~^^ and thus a more severe 
limit on Y~. One finds that the upper limits K™ a * become serious for > 100 GeV and 

lhad Ij ~ 

m f > 0.7 TeV. 

Constraints from late electromagnetic energy injection, Y~^ = £^ x /e em , can become 
more severe than the ones from hadronic energy injection for Tf > 10 4 s. We therefore 
combine the hadronic and electromagnetic constraints to determine the nucleosynthesis 
bounds on the yield of right-handed stau NLSPs before their decay: 

l-BN = min(Yi-d, Y?£) • (4-9) 

In Fig. 12 the resulting severe (left) and conservative (right) bounds — obtained with £h a d 
computed from (3.5) and e em = 0.3 E T (top), 0.5 E T (middle), and E T (bottom) — are shown, 
where the contour lines indicate the upper limits Y-^^ = 10~ 14 , 10~ 13 , 10~ 12 , 10~ n , and 
10~ 10 (as labeled) and a darker shading implies a more severe limit on Yf. The sensitivity 
on e cm is most pronounced for decays at 10 7 s < t < 10 12 s, which are the cosmic times at 
which the constraints from electromagnetic energy release become important. Recall that 
the (electromagnetic) constraints obtained for e em (r — > r G) = E T are overly restrictive for 
stau NLSP scenarios as discussed in Sec. 2. Thus, a stau NLSP scenario which respects 
these constraints preserves the successful predictions of primordial nucleosynthesis. 

Let us comment on the constraints in the alternative selectron NLSP and smuon NLSP 
scenarios. The bounds from hadronic energy injection shown in Fig. 11 can be applied 
directly to eR or /Ir NLSP scenarios by performing the obvious substitutions Yf — > Yg^ and 
mf — > m%,fi- To take into account also the electromagnetic energy injection, we recall from 
Sec. 2 that the electromagnetic energy release from the decay of a single selectron or smuon 
NLSP is described respectively by e em (e eG) ~ £ e and e em (p — ► fxG) ~ 0.3 E^ for > 
1.5 x 10 4 s. Thus, the upper limits on the abundance of selectron and smuon NLSPs before 
their decay — derived from both hadronic and electromagnetic nucleosynthesis constraints — 
can be read respectively from the plots at the bottom and the top in Fig. 12 once the obvious 
substitutions have been performed. 

5. Gravitino Dark Matter from Charged Slepton NLSP Decays 

In this section we consider the non-thermal production (NTP) of gravitino dark matter 
from late decays of charged slepton NLSPs. Additional upper limits on the yield of the 
charged slepton NLSPs before their decay can be obtained since the resulting gravitino 
density cannot exceed the observed dark matter density. We also provide estimates of the 
present free-streaming velocity and the comoving free-streaming scale of gravitino dark 
matter from charged slepton NLSP decays. Associated limits from observations of small- 
scale structure and early reionization are discussed. 

5.1 Relic Gravitino Density from NLSP Decays 

Since each NLSP decays into one gravitino LSP, the resulting relic density of gravitino 
dark matter is governed by the abundance of the NLSPs before their decay. Using h = 
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Figure 12: Upper limits on the yield of right-handed stau NLSPs before their decay, Yf^BN = 
min(l~^^, Y~™£), as obtained from the severe (left) and the conservative (right) upper limits on 
late hadronic and electromagnetic energy injection specified at the end of Sec. 4.1. We computed 
y ?Tad from ( 3 - 5 ) and for £ cm = 0.3 E T (top), 0.5 E T (middle), and E T (bottom). Smaller 

values of Y~^£ N are indicated by a darker shading. 



0.711 'q 3 [52, 7] to parametrize the Hubble constant Hq = 100 /ikm/sec/Mpc, the density 
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parameter of gravitinos produced non-thermally in late NLSP decays reads [13, 15] 

nf p h 2 = mg Fnlsp s(T )h 2 / Pc , (5.1) 

where p c /[s(T )h 2 ] = 3.6 x 10 -9 GeV as obtained from the critical density p c /h 2 = 8.1 x 
1(T 47 GeV 4 , the present temperature T = 2.73 K = 2.35 x 1(T 13 GeV, and the number of 
effectively massless degrees of freedom governing the entropy density today g*s(To) = 3.91. 
The observed density of cold dark matter (95% CL) [52, 7] 1 

^cdm^ 2 = O-llSiaoU (5.2) 
thus gives the following upper limit on the yield of the NLSPs before their decay: 

™ _ n^ M h 2 Pc 12 ( n™ M h 2 \ / loo Gev \ 

y NLSPCDM-/^—^ryp- 4.7x10 f j j (5.3) 

where / denotes the fraction of dark matter which is assumed to be provided by NLSP 
decays, 0~ TP = /^cdm- We- also consider values of / less than one since thermal pro- 
duction of gravitinos in the early Universe can contribute significantly to r^DM' ^ ec - ^- 
Note that (5.1) and (5.3) apply not only to charged slepton NLSPs but also to any NLSP 
decaying into the gravitino LSP in SUSY scenarios with conserved R-parity. 

In Fig. 13 we show ^nlspcdm f° r scenarios with / = 1 (grey band) and / = 0.01 and 
0.1 (solid lines, as labeled). For charged slepton NLSP cases, a comparison with the BBN 
bounds — given in Fig. 12 — shows that the bounds from ^cdm become more severe towards 
larger values of mg, particularly for small values of /. 

5.2 Free Streaming of Gravitinos and Constraints from Small-Scale Structure 

The extremely weakly interacting gravitinos can be treated as collisionless dark matter. 
With a non-negligible velocity, they can stream out of overdense regions and into under- 
dense regions before the time of matter-radiation equality. This free streaming leads to 
an erasure of structure on scales below the comoving free-streaming scale Afs [57]. Thus, 
observations and simulations of cosmic structures can provide severe upper limits on Afs 
or, equivalently, on the present free-streaming velocity Vp S of gravitino dark matter [13, 27]. 
On the other hand, numerical simulations of cold (i.e. t> FS ~ 0) dark matter predict too 
much power on small scales (< 1 Mpc), such as an excess in the number of low mass 
halos and over dense halo cores. These problems could be resolved in warm dark matter 
scenarios (see [69, 62, 63] and references therein) including those in which gravitinos with 
a non-negligible velocity dispersion Vp S dominate ^cdm i^&i 59]. 

For gravitino dark matter produced non-thermally in late NLSP decays, the comoving 
free-streaming scale at the time of matter-radiation equality t eq can be estimated as follows 

FS J tnlsp " " V * ^NLSP 1 + Vl + ^(1+W 2 (WTNLSP) J 

' (5-4) 

lr The value given in (5.2) can be updated taking into account the recent three year results of the Wilkinson 
Microwave Anisotropy Probe (WMAP) [56] . However, this would affect our results at most marginally. 
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Figure 13: Upper limits on the yield of the NLSPs before their decay, i^xsp CDM , from the observed 
density of dark matter ^l^ M h 2 as a function of the gravitino mass mg. The grey band indicates 
the values of (m g , Fnlsp) for which Q^ TP h 2 agrees with Q° b D s M /i 2 = 0.113±g;^| (95% CL) [52, 7]. 
The region above the grey band is excluded. The solid lines show the limits obtained for / ^qdm^ 2 
with n^ x M h 2 = 0.129 and / = 0.01 and 0.1 (as labeled). 



with the cosmic scale factor a(t), the gravitino velocity v(t) = \p(t)\/ \/\p(t)\ 2 + m ~ and 

y G 

its present value vq = t>p S , which is assumed to be non-relativistic. The time and redshift 
at matter-radiation equality are given respectively by 



z eq = 3224 



0.135 



and t cq = 4.7 1 x 1(T 6 = 2.0 x 10 1Z s 



U2 . 



(5.5) 



where to = 4.3 x 10 17 s is the age of the Universe in a cosmological model with matter 
density Q m = 0.27, radiation density JIr = 5 x 10 -5 , and cosmological constant Q\ = 0.73. 
In the estimate (5.4), we use again the sudden decay approximation by assuming that 
the gravitinos are injected when the cosmic time t equals the NLSP lifetime tnlsp- For 
tnlsp "C t eq , the estimate (5.4) agrees with the one given in Refs. [60, 61], 



\GNTP 
A FS 



T"NLSP<*o 



L 



cq , v(t) 
o a(t) 



(l + z eq ) 2 In Jl + - 



1 



(1+z, 



+ 



Vq (1 + Z, 



cq) 



(5.6) 

In charged slepton NLSP scenarios, the non-thermally produced gravitinos originate 
mainly from the 2-body decay 1r — > \G discussed in Sec. 2. The contribution of the 4-body 
decay 1r — ► lGqq to f2~ TP is less than 1% in the considered parameter range as can be 
seen in Figs. 4 and 5. Accordingly, the initial momentum of the gravitino from the 2-body 
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decay tr — > rG 



2 2 2 

m~ — m ~ — mt 



governs the gravitino spectrum, which is monochromatic in the sudden decay (tj = Tf) 
approximation. (A treatment of the phase space distribution beyond the sudden decay 
approximation can be found in Ref. [58].) Under the reasonable assumption that the 
gravitino momentum is non-relativistic today, |pg(io)| = m <g v 0i one can derive directly the 
free-streaming velocity of gravitinos today [27]. With the lifetime Tf given by the inverse of 
the decay width (2.2), we obtain the following result for the present free-streaming velocity 
of gravitinos from stau NLSP decays in the limit m T — > 

<GNTP = Q024 km ( 9 jM=r^y 4 /_m^N-3/2 / _ rr%\- ' 
FS s V 3.36 J VI TeV/ ^ m% ) 1 

where g*(t) is the effective number of relativistic degrees of freedom governing the energy 
density. The result (5.8) is valid for gravitinos from stau NLSP decays in the radiation- 
dominated epoch, Tf < i eq . 

For warm dark matter with a thermal spectrum — such as light (ttiq < 100 eV) grav- 
itinos once in thermal equilibrium [8] — constraints on the associated free-streaming length 
Apg DM and the present free-streaming velocity Upg VDM have been derived from observations 
and simulations of cosmic structures. As demonstrated in [60], these limits can be adopted 
to some extent to non-thermally produced dark matter with a monochromatic spectrum. 
Thus, we summarize the limits obtained for warm dark matter with a thermal spectrum 
and then discuss the implications for gravitino dark matter from charged slepton NLSP 
decays. A similar study with emphasis on constraints from early reionization can be found 
in [27]. 

In warm dark matter investigations of cosmic structure formation, it is typically as- 
sumed that all dark matter consists of a particle species with mass mwDM that was once in 
thermal equilibrium with the primordial plasma and freezes out with a thermal spectrum 
while relativistic at Tf 3> mwDM- For a Majorana fermion of spin 1/2, i.e., a fermion with 
two internal degrees of freedom, the resulting relic density reads 

so that the requirement ^cdm = ^wdm fixes g*s (Tf) for a given w-wdm- This, in turn, 
determines the root mean squared value of the warm dark matter velocity dispersion today 

, rms ,(K n 77 km f n WDM h 2 \ 1/3 / 100eV\ 4/3 

(v FS wdm = 0.77 — — — . (5.10) 

s V °- 113 / V m WDM / 

The obtained upper limits on (wp™ s '°)wDM are thus often expressed as lower limits on 
mwDM' Note that the constraint ^cdm = ^wdm implies g*s(Tf) > 500 for mwDM ^ 
0.5 keV, which is more than about five (two) times the number of degrees of freedom in the 
standard model (MSSM) and requires physics beyond the MSSM (see also the discussion 
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Table 1: Constraints on warm dark matter models from observations and simulations of cosmic 
structures. The lower limits on the mass, m^^ )Ml are taken from the corresponding references. The 
values for (%™ S '°)^J M and (Afs)wdm are derived from to^dm usm g (5.10) and (5.6), respectively. 



™niin / rms,U\ max / \„ n \max 

m WDM \ U FS /WDM l A FSJwDM 

Probe/Observable [keV] [km/s] [Mpc] Ref. 

Dwarf spheroidal galaxies 0.7 0.06 0.64 [63] 

Lyman-Q forest at z ~ 3 0.75 0.05 0.59 [64] 

Lyman-a forest at z ~ (2 - 3) 0.55 0.08 0.84 [65] 

Supermassive black hole at z ~ 5.8 0.5 0.09 0.94 [66] 

Cosmological reionization by z ~ 5.8 0.75 0.05 0.59 [66] 



in [62]). Since if <C t eq , we can use expression (5.6) with vq = (fp™ s ' )wdm to estimate the 
comoving free-streaming scale of the thermal relic at matter-radiation equality, Apg DM . 

In Table 1 we list constraints on w-wdm obtained in warm dark matter studies of 
cosmic structure formation. The considered probes include the maximum observed phase- 
space density in dwarf spheroidal galaxies [63], the observed properties of the Lyman-a 
forest in quasar spectra [64, 65], the supermassive black hole believed to power the quasar 
SDSS 1044-1215, and the reionization inferred from ionized hydrogen in the intergalactic 
medium [66]. With the lower bound on the mass, m^'g M , taken from the given reference, 
the values for (%™ s ' )wdm an d (^fs)wdm were derived from (5.10) and (5.6), respectively. 
The limits from the various probes are basically consistent. One finds that present free- 
streaming velocities of (vpg S '°)wDM ^ (0.05 — 0.09) km/s are disfavored if all dark matter 
consists of warm dark matter with such a velocity dispersion. From cosmological reioniza- 
tion by redshift z = 5.8, a somewhat more severe limit of ttt-wdm ^ 1-2 keV has also been 
discussed in [66], which corresponds to (fp™ s '°)wDM ^5 0.03 km/s. With evidence for cos- 
mological reionization at higher redshift, the limit on (v^ s ' )wdm can become even more 
severe since warm dark matter with a sizeable velocity dispersion delays the formation of 
an early generation of stars or quasars, which provides the most natural explanation for 
early reionization. Indeed, the polarization analysis of the three year WMAP data con- 
strains the optical depth to r = 0.088_ ' 034 , which points to reionization at z = 10.9I 2 ' 3 
in a model with instantaneous reionization [67]. These new results replace the ones from 
the analysis of the first year WMAP data, in which r = 0.17 ± 0.06 was interpreted to 
originate from reionization at z ~ 20 [52]. In the light of the polarization analysis of the 
three year WMAP observations [67], we do not list the (possibly overly restrictive) limit 
mwDM ^ 10 keV, corresponding to (f F ™ s '°)wDM ^5 0.002 km/s, which was derived from the 
first year WMAP observations [68] and discussed in [27] for gravitinos from NLSP decays. 

In Fig. 14 we show the contours of the present free-streaming velocity of gravitinos from 
stau NLSP decays obtained from (5.8) with g*(t > 10 s) = 3.36. The values t>£<f NTP = 
0.002, 0.004, 0.01, 0.04, 0.09, 0.3, and 0.8 km/s (from top to bottom) correspond to the 
present values of the velocity dispersion of warm dark matter with to-wdm = 10, 5, 2, 1, 
0.5, 0.2, and 0.1 keV, respectively. The contours can also be associated with the following 
values of the comoving free-streaming scale at matter-radiation equality: Ap S NTP = 0.03, 
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Figure 14: Contours of the present free-streaming velocities of gravitinos from stau NLSP decays, 
Wpg 3NTP , in the parameter space spanned by and rrif. The indicated values v FS GNTP = 0.002, 
0.004, 0.01, 0.04, 0.09, 0.3, and 0.8 km/s (from top to bottom) can be associated with the present 
velocity dispersion of warm dark matter made of thermal relics of mass towdm = 10, 5, 2, 1, 0.5, 
0.2, and 0.1 keV, respectively. 

0.06, 0.2, 0.4, 1.0, 2.7, and 5.8 Mpc (from top to bottom). 

In scenarios in which basically all dark matter consists of gravitinos from NLSP decays 
(/ ~ 1) — also known as super WIMP gravitino dark matter scenarios [17, 18, 21] — the warm 
dark matter limits from Table 1 can be applied [60]. Accordingly, Fig. 14 shows that such 
scenarios with rrif < 0.5 TeV are in conflict with the constraints inferred from warm dark 
matter investigations of cosmological structure formation. This is a conservative limit. 
With a better understanding of cosmological reionization and the Next Generation Space 
Telescope [66], one could find that even values of the stau NLSP mass up to about 1 TeV are 
excluded. (The additional BBN constraints for super WIMP gravitino dark matter scenarios 
will be updated in Sec. 6.) However, if regions with 0.01 km/s < t> FS GNTP < 0.05 km/s 
remain allowed, superWIMP gravitino dark matter could behave like warm dark matter 
made of thermal relics with 750 keV < w-wdm ^ 2 keV and thus resolve the small-scale 
structure problems of cold dark matter [13, 58, 59]. 

For superWIMP gravitino dark matter from charged NLSP decays at late times, 
TNLSP > 10 7 s, one might think that there is an additional suppression of power on sub- 
galactic scales due to the coupling of the charged NLSP to the photon-baryon fluid [70]. 
However, we will show explicitly that the BBN constraints exclude lifetimes of tnlsp ^ 10 7 s 
completely for the considered charged slepton NLSP scenarios with / ~ 1. 
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For scenarios in which only a fraction / < 1 of fi^DM comprises of gravitinos from 
NLSP decays, the warm dark matter limits listed in Table 1 will be overly restrictive. In the 
limiting case in which the gravitinos behave as hot dark matter, existing constraints from 
studies of hot dark matter such as active neutrinos could be adopted [27]. However, these 
constraints apply only to the high velocity region below the contour ■u F 'g ?NTP = 0.8 km/s in 
Fig. 14, which is excluded anyhow either by the mass bounds from searches for long-lived 
charged particles at LEP [55] or — as will be shown below — by the BBN constraints. In the 
parameter region with ■Upg ?NTP < 0.8 km/s, dedicated N-body simulations are needed to 
extract reliable constraints for such mixed dark matter scenarios from cosmic structures 
on small scales. We postpone this non-trivial task for future work. 



6. Bounds on the Gravitino and Charged Slepton NLSP Masses 

In this section we use estimates of the abundance of charged slepton NLSPs before their 
decay to extract mass bounds for the gravitino LSP and the charged slepton NLSP from 
the cosmological constraints discussed above. Based on our exact treatment of the hadronic 
energy release, we present an update of the mass bounds given in Refs. [15, 21, 22]. We 
consider estimates of the charged slepton NLSP abundance prior to decay, which assume 
that the NLSPs were once in thermal equilibrium and decoupled before decaying into the 
gravitino LSP. In addition, we provide mass bounds for scenarios in which a fixed fraction 
/ of the dark matter density consists of gravitinos from late decays of charged slepton 
NLSPs. 

6.1 Estimates of the Abundance of Charged Slepton NLSPs before their Decay 

The charged slepton NLSPs were in thermal equilibrium with the primordial plasma un- 
til the temperature of the Universe dropped below their freeze-out temperature of Tf ~ 
mj/29 — mj/25 [15]. Since Tf <C mj, these NLSPs are already highly non-relativistic when 
they decouple from the thermal plasma and thus behave as cold thermal relics before their 
decay. This also demonstrates that the charged slepton NLSPs can indeed be considered 
to be at rest when decays into the gravitino LSP at temperatures T <C Tf are studied. 

Although we do not have unambiguous evidence for temperatures of the Universe 
higher than T w 1 MeV (i.e. the temperature at primordial nucleosynthesis), we assume a 
cosmological scenario with temperatures above Tf at early times. This means temperatures 
higher than about at least 3 GeV which is the lowest possible value of Tf obtained from 
the smallest allowed mass of a long-lived charged slepton, my ps 100 GeV [55]. Under this 
assumption, the NLSP abundance becomes independent of the reheating temperature after 
inflation. 

Once the charged slepton NLSPs are decoupled from the thermal plasma, their yield 
lj = nj/s prior to their decay is given approximately by the value at freeze out: Y^(T) ps 
Yp uil (T { ) for T < T f ; cf. [57, 71, 15] and references therein. The precise value of Yy(i = tj) 
depends on details of the NLSP decoupling, which are governed by the mass spectrum 
and the couplings of the superparticles. Indeed, dedicated computer programs such as 
DarkSUSY [72] or micrOMEGAs [73] are available for the computation of the freeze out 



- 26 - 



and the resulting relic abundances of neutralinos and other weakly interacting massive 
particles (WIMPs) in a given SUSY model. 

In this paper we consider the following estimates of the stau NLSP abundance 



which are applicable for nif > 100 GeV as can be seen from the curves in Fig. 1 of Ref. [15]. 
These curves have been derived for = 1.1 m,f and purely right-handed lighter sleptons 
h = 1r (i-e. negligible left-right mixing of sleptons). Coannihilation processes of sleptons 
with binos have not been taken into account. 

The estimate (6.1) is valid for a superparticle spectrum in which the stau NLSP mass 
is significantly below the masses of the lighter selectron and the lighter smuon. Because of 
mfj <C mgj ^ , the abundances of the lighter selectron and the lighter smuon are negligible 
at the time of the stau NLSP freeze out so that selectron/smuon-stau coannihilation pro- 
cesses can be ignored. Accordingly, the selectrons/smuons do not affect the freeze out of 
the stau NLSPs and the stau NLSP abundance is given by Yf ~ rif / s with rif = + rif * . 

The estimate (6.2) is valid for a superparticle spectrum in which the stau NLSP mass 
and the masses of the lighter selectron and smuon are nearly degenerate. Because of 
mfj pa Tn% lt p, 1 , the abundances of the lighter selectron and the lighter smuon are com- 
parable to the stau NLSP abundance and selectron/smuon-stau coannihilation processes 
are important [71, 15]. Although these coannihilation processes reduce the abundance of 
each lighter charged slepton species, the net abundance of the stau NLSP is doubled since 
Yf(rf) ~ Y\(T{) = (ngj + rig* + + n^* + rif x + rif*)/s as each lighter selectron/smuon 
decays eventually into one stau NLSP via ei(/Ii) — > rir/x(e). Note however that these 
decays will be forbidden kinematically if the mass difference — is smaller than 

the mass of the tau lepton. The lighter selectrons/smuons then decay directly into grav- 
itinos and standard model particles. Thus, cosmological constraints have to be considered 
for each lighter charged slepton species separately. This affects the electromagnetic BBN 
constraints as can be seen from our discussion of the selectron/smuon NLSP scenarios in 
Sec. 2. The resulting mass bounds however will remain between the limiting curves shown 
below. In contrast to the electromagnetic BBN constraints, the hadronic ones from late 
( T i ^ 100 s ) decays of charged sleptons are equally valid for each charged slepton species 
(li = 1r) as already mentioned in Sec. 3. The resulting mass bounds will therefore not be 
affected by the discussed precise slepton mass degeneracies. 

In Fig. 15 we plot If as a function of m,f. The solid and dashed curves show the 
estimates (6.1) and (6.2) respectively. Note that the estimate (6.1) has also been used in 
Refs. [22, 74]. It is in good agreement with the result from micrOMEGAs obtained for 
tan/3 = 30 in Ref. [19]. Also the other estimate (6.2) seems reasonable. For example, an 
only slightly higher stau NLSP abundance was found in Ref. [71] with slepton coannihila- 
tions taken into account. 

We consider the estimates (6.1) and (6.2) as representative abundances. The actual 
stau NLSP abundance may be smaller or larger depending on the SUSY model realized 
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Figure 15: The yield of the stau NLSP prior to its decay, Yf = riffs with rif = nf ± + n?*, as 
a function of rrif. The solid and dashed curves show respectively the estimate for <C m^^^ 
given in (6.1) and the one for w m^,^ given in (6.2). Both estimates are taken from Fig. 1 of 
Ref. [15] and have been derived for nig = 1.1 rrif and purely right-handed lighter sleptons, b = 1r. 



in nature. For example, for rrif < < 1.1 mf, the stau-bino coannihilation processes 
become important, thereby leading to an enhancement of the stau NLSP abundance [71]. 
On the other hand, a sizeable left-right mixing of the stau NLSP is associated with an 
increase of its MSSM couplings and thus with a reduction of the stau NLSP abundance [57]. 
Once the SUSY model realized in nature is probed at future colliders, we will be able to 
reduce such uncertainties in the computation of the stau NLSP yield. For a mainly right- 
handed stau NLSP, the resulting yield can then be confronted directly with the cosmological 
upper limits given in Figs. 11-13 above. 

6.2 Mass Bounds from the Estimates of the Charged Slepton NLSP Abundance 

Let us now confront the stau NLSP abundance (6.1) obtained for <C rn^,^ with the 
BBN constraints shown in Fig. 12 and the ^cdm constraints shown in Fig. 13. This leads 
to the bounds on mg and rrif shown in Fig. 16. The observed dark matter density f2^}g M /i 2 
excludes the region above the grey band. The grey band marks the region of superWIMP 
gravitino dark matter scenarios [17, 18, 21] in which fi^DM is provided by gravitinos from 
stau NLSP decays alone, Og TP « ngg M . Only 10% (1%) of f2gg M is provided by gravitinos 
from stau NLSP decays for scenarios that fall onto the thin solid line labeled by / = 0.1 
(/ = 0.01). The thick solid (red) and thick dashed (blue) curves show the limits from late 
hadronic and electromagnetic energy injection, respectively, for both the severe (left plot) 
and the conservative (right plot) BBN bounds. The region to the left of the corresponding 
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Figure 16: Exclusion limits for and m? from the severe (left plot) and conservative (right 
plot) BBN bounds on electromagnetic (thick dashed lines, blue) and hadronic (thick solid lines, 
red) energy release and from ficDM (§ re y band) obtained with the stau NLSP abundance (6.1). 
Allowed is the region from the grey band downwards, to the left of the thick dashed curve(s), and 
outside of the area enclosed by the thick solid curve. On the grey band and on the thin solid 
lines labeled by / = 0.1 and / = 0.01, gravitinos from stau NLSP decays provide 100%, 10%, and 
1% of ^CDM' respectively. The thin dashed lines show contours of Tf — 1 s, 100 s, 10 4 s, and 
10 7 s (from top to bottom). Our study concentrates on Tf > 100 s. The dotted (red) line labeled 
with "FST '04" shows the outdated hadronic BBN constraint from the Chad estimate of Ref. [21]. 
The thick (red) dots indicate the benchmark scenarios A\-E\ that will be introduced in Sec. 9 (cf. 
Table 2). 

curves remains allowed. For the electromagnetic BBN constraints, we obtain the short- 
dashed and the long-dashed curves with the limiting values e em = E T and e em = 0.3 E T 
respectively. The actual BBN exclusion limit from late electromagnetic energy injection 
will fall between these two curves. The dotted (red) line shows the outdated hadronic BBN 
constraint obtained with the simplified estimate (3.6) from Ref. [21]. The thick (red) dots 
indicate the benchmark scenarios A\—E\ that will be introduced and discussed in Sec. 9. 
The thin dashed lines show contours of the stau NLSP lifetime of Tf = 1 s, 100 s, 10 s, and 
10 s (from top to bottom). Recall that we focus our investigation on late decays Tf > 100 s. 
For shorter lifetimes Tf, additional exclusion limits may arise from decays of the tau lepton 
into mesons since these mesons can trigger proton-neutron interconversion processes [29]. 

From a comparison of our new result for the hadronic BBN constraint labeled "had" 
with the outdated one labeled "FST '04", we find that the previous estimate (3.6) from 
Ref. [21] leads to overly restrictive limits. Accordingly, the hadronic BBN constraints given 
in Refs. [22, 23], which were obtained with the estimate of [21], are overly restrictive. 

Our results on hadronic energy release have consequences for the possible existence of 
superWIMP gravitino dark matter scenarios (/ ~ 1) with charged slepton NLSPs. Because 
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of the previously overestimated hadronic energy release, such scenarios were expected to 
be completely excluded for Tf > 100 s by the severe BBN bounds [22]. Also with the 
conservative BBN bounds and particularly for e em ~ E T , such scenarios were believed to be 
almost completely excluded as can be seen in Fig. 16. With our update of the hadronic BBN 
constraints and the yield Y~ n ei ' n given in (6.1), we find that super WIMP gravitino 
dark matter scenarios can be realized around the point D±: (mg,raf) = (200 GeV, 2.8 TeV) 
within the grey region and above the solid (red) hadronic exclusion limit. Here, however, 
one should keep in mind that hadronic tau decays could lead to additional constraints when 
the stau NLSP lifetime approaches 100 s. If the severe BBN bounds are too aggressive and 
the conservative ones are confirmed, superWIMP scenarios will also become allowed within 
the grey band around the point Ey. (m^m?) = (450 GeV, 1.2 TeV), i.e., below the solid 
(red) exclusion limit and above the dashed (blue) exclusion limit. In both of these allowed 
regions, superWIMP gravitino dark matter behaves as warm dark matter (see Fig. 14) 
and could resolve the small-scale structure problems of cold dark matter as suggested in 
Refs. [58, 59]. 

From a comparison of the BBN constraints shown in the left plot with the ones shown 
in the right plot of Fig. 16, one finds a strong sensitivity of the mass bounds on the 
upper limits had snown m Fig- 9. For example, constraints obtained from the severe 
upper limits on £ cm (almost completely) exclude scenarios with lifetimes Tf > 10 7 s. In 
contrast, even significantly larger lifetimes remain allowed with the conservative upper 
limits on £ em . Since the upper limits on late energy injection C^had depend strongly on 
the observed abundances of the primordial light elements, more precise determinations of 
these abundances will be crucial for future refinements of the exclusion limits. 

Note that additional CMB constraints have been extracted from the high precision 
spectral measurements of the Cosmic Background Explorer Far Infrared Absolute Spec- 
trophotometer (COBE FIRAS) [75]. These measurements found a very precise Planck 
spectrum. Possible deviations from such a spectrum can be described by a Bose-Einstein 
distribution with a chemical potential for which |^| < 9 x 10~ 5 has been found [75]. Since 
late electromagnetic energy release can lead to spectral distortions, this upper limit on \fi\ 
can be translated into an upper limit on 6 em [76] . For late decays of the stau NLSP into the 
gravitino LSP, constraints on and m,f inferred from the analytic approximation of [76] 
were given in Refs. [21, 22]. Using a numerical treatment of the kinetic equations for the 
photon number density, these constraints have recently been updated [74] with the estimate 
Y~ T1 1,n given in (6.1). For e em = E T and e em = 0.3 £V, the resulting CMB constraints 
are similar to the corresponding electromagnetic BBN constraints shown in the right plot 
of Fig. 16. Only towards the region around (mg,mf) = (200 GeV, 0.5 TeV), the CMB 
constraint obtained with e cm = E T becomes slightly more severe than the corresponding 
conservative BBN limit. In particular, the 'gap' in the electromagnetic BBN constraint 
for e cm = E T is excluded by the CMB constraint for e em = E T . The electromagnetic BBN 
constraints obtained from the severe BBN bounds however are much more severe than the 
CMB limits. In conclusion, since the present CMB constraints reduce the allowed region 
of (mg,mf) values only slightly when the conservative BBN bounds are considered, these 
limits will not be discussed further in the remainder of this paper. The CMB limits would 
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Figure 17: Same as Fig. 16 but for the stau NLSP abundance (6.2). The thick (blue) dots A 2 -E 2 
indicate the benchmark scenarios that will be introduced in Sec. 9 (cf. Table 3). 



become more important if an improved limit of \fx\ < 2 x 10 could be obtained in future 
missions such as the Absolute Radiometer for Cosmology, Astrophysics Diffuse Emission 
(ARCADE) or the Diffuse Microwave Emission Survey (DIMES) [74]. 

To illustrate the dependence of the mass bounds on the NLSP abundance, we extract 
now the bounds with the stau NLSP abundance (6.2) obtained for ~ mg^^, which 
is twice as large as the one considered in Fig. 16. The results are shown in Fig. 17. For 
the contours and exclusion limits, we use the same labels and plot styles as in Fig. 16. 
The thick (blue) dots indicate the benchmark scenarios A2-E2 that will be introduced and 
discussed in Sec. 9. We find that the mass bounds in Fig. 17 have the same main features 
as the ones in Fig. 16. However, the region of allowed (m^,,rrif) values becomes smaller for 
the larger stau NLSP yield. Moreover, the allowed regions for super WIMP gravitino dark 
matter scenarios move to smaller values of m,Q and mf , i.e., to the regions around the point 
D 2 : (m d ,m f ) = (120 GeV, 2.2 TeV) and the point E 2 : (m d ,mf) = (280 GeV, 1 TeV). 
The latter of these will only be allowed for the conservative BBN bounds, i.e., if the severe 
BBN bounds turn out to be overly restrictive. 

Note that thermal gravitino production can only contribute a negligible amount to the 
relic gravitino density in a superWIMP gravitino dark matter scenario with f^ TP ss ^cdm- 
Thus, for such a superWIMP gravitino dark matter scenario, the upper limits on the 
reheating temperature become severe; see Sec. 7. However, in most of the allowed region 
of (m^, rrif) values, gravitinos from NLSP decays provide only a small fraction of ^cdm so 
that thermally produced gravitinos can be the dominant component of cold dark matter. 

At future colliders, one might be able to identify the point in the (rag, m.f ) plane 
realized in nature. Such prospects and associated implications for cosmology and particle 
physics will be described in Sec. 8. Already before the advent of new collider data, we find 
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upper bounds on the gravitino mass mg of about (100 — 300) GeV and (400 — 700) GeV 
from the exclusion limits in Figs. 16 and 17 respectively. These bounds can be refined 
and possibly tightened by computing the precise value of e cm and by extending our study 
to smaller stau NLSP lifetimes. This is an important task since an upper limit on mg 
constrains the SUSY breaking scale which is crucial for insights into the SUSY breaking 
mechanism. 



6.3 Mass Bounds for Gravitino Dark Matter Fractions from Slepton NLSP 
Decays 

We now consider the mass bounds in scenarios in which a fixed fraction / of the observed 
amount of dark matter consists of gravitinos from late decays of charged slepton NLSPs. 
Such scenarios require an NLSP yield prior to decay of 

W = / * = 4.1 x 10-" / (%f) (HgV) , (6.3, 

s(T )h z \ 0.113 / \ / 

so that the bound from is respected by definition; see also Fig. 13. Accordingly, 

one is in the allowed region whenever the ^cdm constraints are more severe than the BBN 
constraints shown in Fig. 12. Note that the abundance (6.3) depends only on mg for a 
given /. In contrast, the estimates of the thermal relic abundance given in (6.1) and (6.2) 
depend on mf but are independent of uiq. Indeed, only if the abundance (6.3) agrees with 
the estimates of the thermal relic abundance, will we consider it to be in a natural range. 

In Fig. 18 we show the exclusion limits on mg and nif obtained when (6.3) with / = 1, 
0.1, and 0.01 (from top to bottom) is confronted with the severe (left) and conservative 
(right) BBN bounds for late (rf > 100 s) NLSP decays. Contours of Tf are also shown. 
For the BBN constraints and the Tf contours, we use the labels and plot styles used al- 
ready in Figs. 16 and 17. The allowed region is the one to the left of the curves from 
electromagnetic energy release (thick dashed lines, blue) and outside of the parameter re- 
gion bounded by the curves from hadronic energy release (thick solid lines, red). Recall 
however that we consider scenarios with Tf > 100 s and that additional constraints can 
appear for Tf < 100 s. The grey bands indicate the regions in which the corresponding 
NLSP abundance (6.3) seems natural as it falls between the estimates of the thermal relic 
abundance given in (6.1) and (6.2). The exclusion limits inside the grey bands therefore 
resemble the ones found in Figs. 16 and 17. Above the grey bands, the abundance (6.3) 
becomes smaller than Y™ 7 ^ 1 m ° 1,Al S o that (mg,mf) values are allowed that should be 
excluded according to Fig. 16. Below the grey bands, the abundance (6.3) becomes larger 
than Y™ T1 ~ m,11 ' n so that the BBN constraints exclude (mg,mf) values in Fig. 18 which 
are allowed in Fig. 17. For the considered scenarios, one sees clearly the increase of the 
allowed (m^,mf) region with decreasing /. Moreover, the dotted (red) lines show that 
the hadronic BBN constraints were overestimated towards large values of mf and under- 
estimated towards small values of mf by the previous eh a d estimate of Ref. [21]. This is 
exactly the failure of the simplified treatment expected from Fig. 8. 
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Figure 18: Exclusion limits for and m? from the severe (left) and conservative (right) BBN 
bounds on electromagnetic and hadronic energy release for scenarios in which gravitinos from stau 
NLSP decays provide 100%, 10%, and 1% of n$g M ( from to P to bottom). For the BBN constraints 
and the 7> contours, we use the same labels and plot styles as in Fig. 16. Allowed is the region to 
the left of the thick dashed curve(s) and outside of the area enclosed by the thick solid curve. On 
the grey band, the stau NLSP abundance (6.3) has a natural value bounded by (6.1) and (6.2). 

7. Gravitino Dark Matter from Thermal Production 



In this section we consider the thermal production of gravitinos in the early Universe. We 
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compute the present free-streaming velocity and the comoving free-streaming scale of ther- 
mally produced gravitinos. We present upper limits on the reheating temperature after 
inflation, Tr, for scenarios in which a given fraction of the observed dark matter comprises 
of gravitinos from thermal production. We stress that these limits can be tightened by the 
cosmological constraints discussed in the previous sections. To illustrate the significance 
of the Tr bounds, we address the implications for thermal leptogenesis in a SUSY frame- 
work. A similar discussion can be found in the framework of the CMSSM with a gravitino 
LSP [23]. In our study we are exploring a wide range of the gravitino mass, the slepton 
NLSP mass, and the gluino mass, and we do not restrict our investigation to a constrained 
framework. 

7.1 Relic Gravitino Density from Thermal Production in the Early Universe 

Because of their extremely weak interactions, the temperature Tf at which gravitinos de- 
couple from the thermal plasma can be very high depending on the gravitino mass. For 
m G ^ 10 keV (10 MeV), one finds a gravitino decoupling temperature of Tf > 100 TeV 
(10 8 TeV). However, for a reheating temperature above the decoupling temperature, 
Tr > Tf, the relic density of a gravitino LSP with 10 keV < < 1 TeV exceeds 
^cdm significantly for any reasonable value of g*s(T{) such as g*s (Tf) — 230, as can 
be seen from (5.9). Thus, only gravitino LSP scenarios with Tr <C T are allowed for 
10 keV <mg<l TeV. 

Note that a light gravitino LSP with ~ 100 eV has a much smaller decoupling tem- 
perature of Tf = 0(1 TeV). For Tr > Tf, such a gravitino — once in thermal equilibrium — 
would even provide the right amount of dark matter as can be seen from (5.9) with 
9*s(Tf) ~ 100. However, the velocity dispersion of such a light gravitino given by (5.10) 
exceeds significantly the constraints from observations and simulations of cosmic structures 
listed in Table 1. 

For Tr <C Tf, gravitinos can be produced in thermal reactions in the hot MSSM 
plasma. The thermal gravitino production rate at high temperatures is computed in a 
gauge-invariant way [16] with the Braaten-Yuan prescription [77] and hard thermal loop 
resummation [78] , which takes into account Debye screening in the plasma. Assuming that 
inflation has diluted away any primordial gravitino abundance, gravitino disappearance 
processes are negligible for Tr <C Tf. The corresponding Boltzmann equation can be 
solved analytically. To leading order in the gauge coupling of quantum chromodynamics, 
this leads to the following result for the relic density of stable LSP gravitinos from thermal 
production (TP) [16] 

O W - 12 o 2 In (h™) (l + ™L] ( m Q \ ( — ^ (7 1) 

G n g j y- + 3m 2 s J lioOGeVv 1 \l0™ GeV J ' U > 

where g is the strong coupling and m g the gluino mass. The contribution proportional 
to mj results from the spin-1/2 components of the gravitino. Indeed, this contribution 
governs 0^ p /i 2 for mg <C nig. For larger values of the gravitino mass, mg > 0.5 m g , the 
contribution of the spin-3/2 components — which is independent of m g — becomes equally 
important. 
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Figure 19: The relic density of gravitino LSPs from thermal production, Qg P /i 2 , as a function of 
T R for m s (M z ) = 0.7 TeV (top), 1.4 TcV (middle), and 5 TeV (bottom). The dashed lines show the 
results for m g = 10 MeV, 100 MeV, 1 GeV, 10 GeV, 100 GeV, and 400 GeV (from the upper left to 
the lower right). The grey band indicates the observed dark matter density Q(^ M h 2 = O.llS+oloi! 
(95% CL) [52, 7]. On the solid lines labeled with (1 - /) = 0.1 and (1 - /) = 0.01, the contribution 
of thermally produced gravitinos to ^cdm i s 10% an( l 1% respectively. 



In Fig. 19 the dashed lines show VL^h 2 as a function of the reheating temperature 
Tr for different values of the gravitino mass ranging from 10 MeV to 400 GeV (from 
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the upper left to the lower right), where the values rrig(Mz) = 0.7 TeV (top), 1.4 TeV 
(middle), and 5 TeV (bottom) are considered. The running of the strong coupling and 
the gluino mass is taken into account by replacing g and m g in (7.1) respectively with 
g{T R ) = [g- 2 (M z )+3ln(T R /M z )/(8Tr 2 )}~ 1 / 2 and m g {T R ) = [g(T R )/g(M z )] 2 m g (M z ), where 
g 2 (M z )/(4Tr) = 0.118. The upper two plots show explicitly that a doubling of the gluino 
mass enhances Vt^h? by a factor of four for mg <C m g . The results obtained for mg = 
400 GeV demonstrate that the contribution of the spin-3/2 components becomes important 
for mg > 0.5 m g . On the grey band, fiT p agrees with ^cdm- On ^ ne son d lines labeled with 
(1 — /) = 0.1 and (1 — /) = 0.01, the thermally produced gravitinos provide respectively 
10% and 1% of ngg M , which leaves room for a significant contribution of gravitinos from 
stau NLSP decays. 

In studies of scenarios with gravitino LSP and stau NLSP, a gluino mass above the 
mass of the lighter stau has to be considered. For such scenarios, one should also keep in 
mind that a gravitino LSP with mg > 10 GeV can be excluded for certain values of the 
stau NLSP mass by the cosmological constraints shown in Figs. 16 and 17. In particular, 
gravitino masses above about 700 GeV seem to be excluded for any value of the stau NLSP 
mass as stressed at the end of Sec. 6.2. 



7.2 Free Streaming of Gravitinos from Thermal Production 

For a reheating temperature below the gravitino decoupling temperature, gravitinos have 
never been in thermal equilibrium with the primordial plasma. However, since the thermal 
production proceeds through reactions of particles in thermal equilibrium, the gravitinos 
are produced in kinetic equilibrium with the primordial plasma. Therefore, the gravitinos 
from thermal production have a thermal spectrum. The root mean squared value of their 
velocity dispersion today is accordingly given by 



rms.chGTP _, Sy1ft - 6 km /10MeV\ / 230 



;r:r-'- < r , <' - o.s x i.r" — ( ) ( ) 1/3 . (7.2) 

FS s \ m d J Ws(T R )J 

Since the thermal gravitino production is efficient only in the very early hot Universe, the 
corresponding comoving free-streaming scale can be estimated from expression (5.6) with 
v o = (wps S '°) GTP . Thermally produced gravitino LSPs with mg > 100 keV have a negligi- 
ble free-streaming behavior and thus can be classified as cold dark matter. Accordingly, in 
scenarios in which basically all dark matter is made of such gravitinos, one must face the 
small-scale (< 1 Mpc) structure problems inherent to cold dark matter models. 



7.3 Upper Limits on the Reheating Temperature 

Since lH p depends on the reheating temperature after inflation, the observed dark matter 
density allows us to extract upper limits on T R as a function of mg in gravitino LSP sce- 
narios with a given gluino mass [11, 15, 16, 23, 25]. These limits will become particularly 
severe if m g 3> mg and if gravitinos from NLSP decays provide already a sizeable fraction 
/ of ^cdm- The range of allowed values of the reheating temperature is crucial for our 
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understanding of inflation. Moreover, thermal leptogenesis provides an attractive expla- 
nation of the baryon asymmetry in the Universe for very high reheating temperatures of 
T R > 3 x 10 9 GeV [47, 48]. 

In Fig. 20 we show the reheating temperatures Tr for which = (1 — /) ^cdm w ^h 
(1 — /) = 1, 0.1, and 0.01 (as labeled) as a function of mg. The results are presented for 
gluino masses of m g (Mi) = 0.7 TeV (top), 1.4 TeV (middle), and 5 TeV (bottom). Since 
our focus is on gravitino LSP scenarios, we consider only gravitino masses that are smaller 
than the running gluino mass. For the (wig, Tr) values within the grey band, SlT p agrees 
with the observed dark matter density ^cdm- The region above the grey band is excluded 
for any value of the NLSP mass. If the NLSP decays provide already 90% or 99% of figg M , 
we can exclude the region above the solid lines labeled respectively with (1 — /) = 0.1 
or (1 — /) = 0.01. Also the bounds presented for charged slepton NLSP scenarios in 
Figs. 16 and 17 can tighten the Tr limits for > 10 GeV, even in scenarios with / <C 1. 
We use again thick dots to indicate the benchmark scenarios ^.(1,2)^(1,2) (red, blue) that 
will be introduced and discussed in Sec. 9. 

Figure 20 shows that thermal leptogenesis cannot explain the observed value of the 
baryon asymmetry in a SUSY framework with gravitino LSP and a heavy gluino rrig(Mz) = 
5 TeV. Similar investigations have found that successful thermal leptogenesis (Tr > 
3 x 10 9 GeV) implies upper bounds on the gluino mass of about 1.3 TeV and 1.8 TeV 
respectively for charged slepton NLSP and sneutrino NLSP scenarios [19]. We postpone 
an update of the gluino mass bounds based on our improved cosmological constraints for 
future work. 

8. Prospects for Collider Phenomenology and Cosmology 

In this section we describe the possible interplay between astrophysical investigations and 
studies at future particle accelerators — such as the Large Hadron Collider (LHC) and the 
International Linear Collider (ILC) — for gravitino LSP scenarios with a charged slepton 
NLSP. We review potential collider signatures of gravitino dark matter and stress the 
implications for cosmology. 

For phenomenology at future colliders, the SUSY scenario with a charged slepton NLSP 
and an extremely weakly interacting LSP such as the gravitino is particularly promis- 
ing [36, 37, 38, 39, 40, 41, 42, 22, 43, 44, 45, 25]. Because of its long lifetime, the charged 
slepton NLSP would appear as a (quasi-)stable particle. Such a long-lived heavy charged 
particle would penetrate collider detectors in a way similar to muons [36, 37, 39]. A sig- 
nificant fraction of the NLSP decays will take place outside the detector and will thus 
escape detection. If the produced charged slepton NLSPs are slow, the associated highly 
ionizing tracks and time-of-flight measurements will allow one to distinguish the NLSP 
sleptons from muons [36, 37, 39, 41]. From measurements of the velocity /3j = v~Jc and 
the momentum = of the long-lived slepton, its mass can be determined from the 
standard relation mj = P[^J^ — fif/P\ [41]- If some of the charged slepton NLSPs would 
decay already in the detectors, the statistical method proposed in Ref. [41] could allow one 
to measure the lifetime of the slepton NLSP. Since this lifetime is governed by the 2-body 
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Upper Limits on the Reheating Temperature 
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Figure 20: The reheating temperatures Tr for which thermally produced gravitino LSPs provide 
100% (grey band), 10% (middle solid line), and 1% (lower solid line) of figg M ft, 2 = 0.1131^1 
(95% CL) [52, 7], where m s (M z ) = 0.7 TeV (top), 1.4 TeV (middle), and 5 TeV (bottom). The 
region above the grey band is excluded for any value of the NLSP mass. For a given superparticle 
spectrum, an additional upper bound on mg, is obtained which can lead to a more severe upper 
limit on Tr] see Figs. 16 and 17 for stau NLSP scenarios. The thick dots indicate the benchmark 
scenarios Am j) — -^(1,2) ( re d, blue) that will be introduced in Sec. 9 (cf. Tables 2 and 3). 

decay 1 — > G\ in the gravitino LSP case, we can use the decay rate (2.2) with the measured 
value of mj to determine also the gravitino mass mg. This would be an experimental 
determination of the (rog, mr) point in Figs. 16 and 17. 



- 38 - 



In addition, there have been works, which propose ways to stop and collect charged 
long-lived particles for an analysis of their decays [79, 43, 44]. It was found that up to 
0(1O 3 -1O 4 ) and 0(1O 3 -1O 5 ) of charged slepton NLSPs can be trapped per year at the 
LHC and the ILC, respectively, by placing 1-10 kt of massive additional material around 
planned collider detectors [43, 44]. With O(10 4 ) of analyzed NLSP decays, one could 
distinguish between the gravitino LSP scenario and the alternative axino LSP scenario by 
considering the differential distributions of the visible decay products in the 3-body decay 
of the slepton NLSP into the LSP, the corresponding lepton, and a photon [45, 25]. The 
axino is another extremely weakly interacting LSP candidate, which can also be produced 
in decays of charged slepton NLSPs [80, 45, 25] and in thermal reactions in the very early 
Universe [81, 82]. In the axino LSP case, the analysis of the slepton NLSP decays could 
allow one to determine the axino mass and to estimate the Peccei-Quinn scale [45, 25]. 
In the gravitino LSP case, could be inferred kinematically. If so, one will be able 
to use the decay rate (2.1) to determine the Planck scale Mpi microscopically from the 
measured slepton NLSP lifetime and the measured values of mg and mj. If consistent with 
macroscopic measurements of Mpi, this will be strong evidence for the existence of SUGRA 
in nature [42]. 

The prospects for phenomenology with long-lived charged NLSP sleptons at future 
colliders depend strongly on the mass spectrum of the superpartners. The heavier the 
superparticles the smaller will be the total SUSY production cross section and, in particular, 
the charged slepton NLSP production cross section. If mj is smaller than about 0.2 TeV, 
the charged slepton NLSP will be copiously produced in the proton-proton collisions of 
cm. energy y/s = 14 TeV at the LHC and even in the electron-positron collisions of cm. 
energy \fs = 0.5 TeV at the ILC. A long-lived charged slepton NLSP with mj = 0.7 TeV 
could possibly still be discovered at the LHC with a luminosity of 100 fb _1 [22]. At the 
ILC this will require the second phase in which a cm. energy of \fs = 1 TeV might be 
reached. Heavier NLSPs will be difficult to explore at the LHC since the production rate 
becomes very small. Even smaller production rates are expected at the ILC; cf. [22, 43, 44] 
and references therein. 

The appearance of a long-lived charged slepton as the lightest standard model super- 
partner in the collider detector would already be a strong hint towards physics beyond the 
MSSM. Indeed, from the severe constraints for stable charged massive particles (cf. [83] 
and references therein), one would expect that the charged slepton will eventually decay. 
A long-lived charged slepton thus points to an extremely weakly interacting LSP such as 
the gravitino or the axino. 

Once the mass of the long-lived slepton is measured, the cosmological constraints 
shown in Figs. 16 and 17 at the measured value of my could tighten the upper bound 
on the gravitino mass discussed at the end of Sec. 6.2. This upper bound on mg, which 
implicitly assumes the gravitino LSP scenario, will imply an upper bound for the SUSY 
breaking scale. In addition, this constraint will give an upper limit on the amount 
of dark matter that could come from late decays of the charged NLSP sleptons. If it 
turns out that this amount is only a minor fraction (/ < 0.1) of the observed dark matter 
density, there will be room for a significant contribution of thermally produced gravitino 
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dark matter and, thus, for high values of the reheating temperature. 

If one succeeds to analyze the decays of the long-lived sleptons and to find evidence for 
the gravitino LSP, this will provide an experimental determination of the (mg,mj) point 
as described above. The obtained values of mg and will then fix the maximum amount 
of gravitino dark matter from slepton NLSP decays. Since the remaining fraction of the 
dark matter density can be made of thermally produced gravitinos, this will give an upper 
limit on the reheating temperature after inflation, which is crucial for our understanding of 
inflation and the baryon asymmetry of the Universe. At this point, it will also be interesting 
to see whether the obtained (mg, mj) point falls into the region which we expect to be 
allowed by the cosmological constraints shown in Figs. 16 and 17. 

We expect that the cosmological constraints can be refined with future astrophysical 
investigations and insights from future colliders. By reducing the uncertainties in the 
observed abundances of the light primordial elements and by improving the numerical 
studies of the effects of late decaying particles on those elements, one should be able to 
reduce the uncertainty between the limits from the severe and the conservative BBN bounds 
shown in Figs. 16 and 17. If the masses of the lightest SUSY particles and their couplings 
can be determined at future colliders, one will also be able to compute the yield of the 
slepton NLSP before its decay with high precision, particularly, for the scenario in which 
the slepton NLSPs were once in thermal equilibrium. Such a precise determination of Inlsp 
will reduce the uncertainty between the constraints shown in Fig. 16 and the ones shown 
in Fig. 17. Experimental insights into the mass hierarchy of the lightest superpartners can 
also be necessary in order to calculate reliably the electromagnetic energy release e em and 
the associated exclusion limits. 



9. Benchmark Scenarios with Gravitino LSP and Stau NLSP 

To summarize the insights gained from the cosmological constraints presented above, we 
consider ten benchmark scenarios A\^ to £1,2 for the case of the gravitino LSP and the stau 
NLSP with mg not smaller than 10 GeV and = 1.1 mf. We present the scenarios A\ 
to E\ for mfj <C rrie l! p, 1 in Table 2 and the scenarios A2 to E2 for ~ in Table 3. 

Each scenario is given by (m^,mf) with <C or « ^ei,Ai and represented 

respectively by a thick (red) dot in Fig. 16 or by a thick (blue) dot in Fig. 17. We also 
use thick (red, blue) dots to indicate the benchmark points A\^.-E\^ in Fig. 20. In the 
tables, we list for each scenario the lifetime of the stau NLSP Tf , the fraction / of f^DM 
which is provided by gravitinos from stau NLSP decays, the present velocity dispersion 
of such non-thermally produced gravitinos ( - U]t<s) G ' NTP j anc ^ an associated classification 
as a hot, warm, or cold dark matter component. In addition, we indicate whether the 
severe and/or conservative BBN bounds are respected. We assume that the remaining 
fraction (1 — /) of f^DM consists of thermally produced gravitinos and specify also the 
present velocity dispersion of the thermally produced gravitinos (fp S ) GTP which allows us 
to classify this component as cold dark matter. Moreover, the reheating temperature Tr 
required for SlT p = (1 — /)fi^ M is given for gluino masses of rrig(Mz) = 0.7 TeV, 1.4 TeV, 
and 5.0 TeV. Note that higher values of Tr are excluded while lower values of Tr are 
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Table 2: Benchmark scenarios for the gravitino LSP and stau NLSP case with m fl <C m 5l ^ 



Scenario 


A 

A\ 




D\ 


w 


L)\ 


171 

hj\ 


mg [GeV] 


10 




50 


100 


200 


450 


m,f [TeV] 


0.2 




0.2 


0.7 


2.8 


1.2 


Tf [S] 


1.9 x 


10 5 


6.0 x 10 6 


3.8 x 10 4 


140 


8.8 x 10 4 


conservative BB1N bounds 


yes 




yes 


yes 


yes (?) 


yes 


severe BB1N bounds 


yes 






yes 


yes (?) 


no 


/ 


3.5 x 


10~ 3 


1.8 x 10~ 2 


0.12 


0.99 


0.96 


Ks) 6NTP [km/s] 


0.26 




0.28 


0.04 


0.005 


0.02 




hot 




hot 


warm 


warm 


warm 


(C) 6TP [km/s] 


5.8 x 


io- 9 


1.2 x 10~ 9 


5.8 x 10~ 10 


2.9 x 10~ 10 


1.3 x 10~ 10 




cold 




cold 


cold 


cold 


cold 


T ™ 5 =0.7TcV [QeV] 


1.0 x 


10 9 


4.8 x 10 9 


7.3 x 10 9 






T m 5 = 1.4TcV [QeV] 


2.3 x 


10 8 


1.2 x 10 9 


2.0 x 10 9 




2.2 x 10 8 


^=5.0 TeV [QeV] 


1.7 x 


10 7 


8.8 x 10 7 


1.6 x 10 8 


3.0 x 10 6 


3.2 x 10 7 



Table 3: Benchmark scenarios for the gravitino LSP and stau NLSP case with w m^^j. 



Scenario 


A 2 




B 2 


c 2 


D 2 


E 2 


m d [GeV] 


10 




50 


100 


120 


280 


m? [TeV] 


0.2 




0.2 


0.7 


2.2 


1.0 


Tf [s] 


1.9 x 


10 5 


6.0 x 10 6 


3.8 x 10 4 


140 


8.8 x 10 4 


conservative BBN bounds 


yes 




yes 


yes 


yes (?) 


yes 


severe BBN bounds 


yes 




? (Cem) 


yes 


yes (?) 


no 


/ 


7.1 x 


10~ 3 


3.3 x 10~ 2 


0.25 


0.93 


0.99 


« S ) GNTP [km/s] 


0.26 




0.28 


0.04 


0.005 


0.02 




hot 




hot 


warm 


warm 


warm 


(C) 6TP [Ws] 


5.8 x 


io- 9 


1.2 x 10~ 9 


5.8 x 10~ 10 


2.9 x 10~ 10 


1.3 x 10~ 10 




cold 




cold 


cold 


cold 


cold 


T ™ 5 =0.7TcV [QeV] 


1.0 x 


10 9 


4.6 x 10 9 


6.0 x 10 9 






T m g =1.4TeV ^ 


2.3 x 


10 8 


1.2 x 10 9 


1.7 x 10 9 




4.0 x 10 7 


^=5.0 TeV [QeV] 


1.7 x 


10 7 


8.5 x 10 7 


1.3 x 10 8 


1.3 x 10 7 


4.0 x 10 6 



possible with additional gravitino sources such as innaton decays [26] and/or additional 
dark matter constituents such as axions [84]. Depending on the value of /, we can group 
the scenarios to £1,2 into the three classes which are discussed in the following. 

In the scenarios and B\^, the stau NLSP decays contribute less than 5% of the 
observed dark matter density. This minor fraction can be classified as hot dark matter as 
it has a relatively high velocity dispersion today, (wfs) GNTP ^ ^ ^ km/s. The dominant 
contribution to ^qdm can thus cons i s t of thermally produced gravitinos. These gravitinos 
have a negligible velocity dispersion today and thus can be classified as cold dark matter. 
Accordingly, the small-scale structure problems of cold dark matter remain to be resolved. 
Reheating temperatures of Tr > 10 9 GeV will be allowed in these scenarios if the gluino 
mass is not much larger than 1 TeV. In particular, for the points B\p and a gluino mass 
not much larger than 0.7 TeV, thermal leptogenesis provides a viable explanation of the 
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baryon asymmetry. The severe and conservative BBN bounds are respected by the points 
Ai ; 2- The points B\2 respect the conservative BBN bounds as well but will respect the 
severe BBN bounds only if e cm ~ 0.3 E T . For a larger electromagnetic energy release, the 
points B\ % 2 will not be allowed by the severe BBN bounds. The appealing feature of the 
points Ai t 2 and is the relatively low stau NLSP mass of rxif = 0.2 TeV, which would 
make these scenarios accessible at the LHC and also at the ILC already in its first phase 
with \/s = 0.5 TeV. Note that the stau NLSP mass is only four times larger than the 
gravitino mass in the scenarios Bi2- The spin-3/2 components of the gravitino emitted 
in the stau NLSP decay would therefore be non-negligible and even a measurement of the 
gravitino spin would be conceivable [42] . 

In the scenarios Ci ; 2, the gravitinos from stau NLSP decays provide a considerable 
fraction of about (10 — 25)% of the observed dark matter density. This fraction has a non- 
negligible velocity dispersion today, (wfs) GNTP > 0.04 km/s, and can be classified as warm 
dark matter. If the remaining part of ^qdm * s ma de of thermally produced gravitinos, 
which would have a negligible velocity dispersion today, one will have mixed gravitino dark 
matter with a cold and a warm component. Again it is not clear at this point if the small- 
scale structure problems of cold dark matter can be overcome in such a scenario. However, 
such a scenario is clearly allowed by the severe and the conservative BBN constraints. 
At the points Ci 2, the reheating temperature needed for successful thermal leptogenesis, 
T R > 3 x 10 9 GeV, will be allowed if the gluino mass is not much larger than 1 TeV. Note 
the exceptional case of a gluino with a mass of 0.7 TeV which would be as be as light as 
the stau NLSP. In this case, we implicitly assume a finite mass difference since the gluino 
NLSP scenario is excluded by hadronic BBN constraints [19]. A stau NLSP with a mass 
of rrif = 0.7 TeV will be very hard to produce in the LHC experiments and will be out of 
reach in the first phase of the ILC. However, as discussed in Sec. 8, even a small number 
of long-lived stau NLSPs could give an exciting signature and could allow us to determine 
the stau NLSP mass. 

In the scenarios D12 and £1,2, more than 90% of the observed cold dark matter 
density is made of gravitinos from stau NLSP decays. The emitted gravitinos still have 
a non-negligible velocity dispersion today, (^fs) GNTP = 0-005 — 0.02 km/s. As warm 
dark matter, they could explain the observed power on small scales, which is significantly 
lower than expected from simulations of cold dark matter. Since the contribution of the 
thermally produced gravitinos cannot exceed 10% of the observed dark matter density, the 
reheating temperature after inflation is constrained to values below about Tr = 10 s GeV. 
Accordingly, thermal leptogenesis seems excluded as a possible explanation of the baryon 
asymmetry in these scenarios. However, it should be stressed that the scenarios D\ 2 and 
E\ t 2 can be in conflict with the BBN bounds. In particular, the points £1,2 are already 
excluded by the severe BBN bounds. The points D12 seem to be allowed but these points 
are associated with stau NLSP lifetimes close to 100 s, where additional constraints could 
arise from hadronic decays of tau leptons. Unfortunately, the scenarios D\ 2 and £1,2 will 
be extremely hard to probe at the LHC and impossible to probe at the ILC. A stau NLSP 
heavier than about 1 TeV implies that all standard model superpartners are heavier than 
about 1 TeV and, thus, a very small value of even the total SUSY production cross section. 
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Finally, we stress that scenarios with mg < 10 GeV and mf < 5 TeV cannot be 
excluded by the constraints from BBN and ^cdm- ^ n this region of the parameter space, 
the contribution of gravitino dark matter from NLSP decays is below 10% of ^cdm so that 
almost all of the observed amount of dark matter can be explained by thermally produced 
gravitinos. The upper bounds on the reheating temperature however become more severe 
towards smaller values of mg as shown in Fig. 20. Moreover, for mg < 100 keV, the present 
velocity dispersion of thermally produced gravitinos becomes non-negligible so that these 
gravitinos can be classified as warm dark matter. The warm dark matter constraints listed 
in Table 1 can then be used to exclude thermally produced gravitinos with mg < 100 eV. 

10. Conclusions 

We have presented a systematic and thorough investigation of cosmological constraints 
in i?-parity respecting SUSY scenarios with a gravitino LSP and a long-lived (tj > 100 s) 
charged slepton NLSP. The gravitino appears naturally as the LSP in well-motivated SUSY 
breaking schemes [1, 3, 5] and the lighter stau as the lightest standard model superpartner 
in high-energy frameworks. Scenarios with the gravitino LSP and the stau NLSP are 
particularly promising: The gravitino LSP is a natural dark matter candidate and the long- 
lived stau NLSP could provide striking signatures at future colliders [36, 37, 38, 39, 40, 41] 
which could lead to experimental evidence for gravitino dark matter [42, 45, 25]. 

In our study we have assumed that gravitino dark matter originates from thermal 
production in the early Universe [11, 12, 14, 16] and from late NLSP decays [13, 15, 17, 
18, 21, 22], both of which are guaranteed production mechanisms. We have examined the 
relic density of gravitinos, Qq, their present free-streaming velocity, fpg , the release of 
electromagnetic and hadronic energy in late decays of charged NLSP sleptons, e emj h a d> an d 
the yield of the NLSP sleptons before their decay, lj. These quantities are constrained by 
the observed dark matter density, by observations and simulations of cosmic structures, 
and by the observed abundances of the light primordial elements. We have used these 
cosmological constraints to extract bounds on the gravitino mass, mg, the mass of the 
charged slepton NLSP, my, and the reheating temperature after inflation, Tr. In addition, 
upper limits on the yield of purely right-handed NLSP sleptons before their decay, lC max ; 
have been obtained. 

New results for the hadronic energy release ehad in late (tj > 100 s) decays of charged 
right-handed slepton NLSPs have been presented from our computation of the 4-body 
decay of the slepton NLSP into the gravitino, the corresponding lepton, and a quark- 
antiquark pair. In our computation, we have taken into account the exchange of virtual 
electroweak gauge bosons and the exact energy spectrum of the emitted quark-antiquark 
pair. In the previous estimate, only Z-boson exchange in the zero-width approximation was 
considered together with a simplified monochromatic approximation of the energy of the 
quark-antiquark pair [21]. From a comparison with our exact treatment, we find that ehad 
has been underestimated towards smaller values of rrif and overestimated for wif 3> mg,. 
Therefore, the hadronic nucleosynthesis constraints given in Refs. [21, 22, 23], which were 
derived from the simplified estimate of Ref. [21], have to be updated. 
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With our exact computation of e^ad) we have derived new upper limits on from the 
recent nucleosynthesis constraints on late energy injection, £ e m,had(= eem,had ^) , given in 
Ref. [29, 31]. We have applied the BBN bounds from the observed abundances of primordial 
D and 4 He since they are the most reliable ones. To account for the present uncertainties 
of these abundances, we have used one severe and one conservative BBN bound for both 
hadronic and electromagnetic energy injection. We find that the hadronic limits ^j™ 3 ^ 
become serious for mg > 100 GeV and nif > 0.7 TeV, and already for NLSP lifetimes of 
100 sec < tj < 10 4 sec. In contrast, the electromagnetic limits Y- 1 ^^ can be important for 
smaller values of mg and nif but only for NLSP lifetimes of tj > 10 4 sec. 

We have also discussed the upper limit on Yy from the observed dark matter density, 
^cdm> w hi cn arises since the gravitino density from NLSP decays, ^^ TP > cannot exceed 
^cdm- This upper limit ^y™^ M is independent of nif. It becomes severe towards larger 
values of mg and for scenarios in which the thermally produced gravitino density f2^ TP 
dominates ^cdm- 

The upper limits Y- max , K™ a *, and K™^,, have allowed us to extract bounds on 

1 ^ lcm ' lhad' 1CDM Ct 

and nif for scenarios in which the NLSP sleptons freeze out with a thermal abundance. 
Since the thermal relic abundance of the slepton NLSP prior to decay is sensitive to the 
mass spectrum and the couplings of the superparticles, we have considered an estimate 
Y™ ri mci,iil for mfj <C and another estimate Y™ T1 moi,Al for ~ mg^/^, 

where Y™ T1 mci,iil [ s two times larger than Y™ T1 m ° 1 '' 31 because of selectron/smuon-stau 
coannihilation processes [15]. In both estimates, the yield Yf is proportional to nif. We 
have indicated typical uncertainties of the mass bounds with respect to the SUSY model by 
working with both Yf estimates. Another uncertainty is related to the precise value of 6 em 
which can also depend on the superparticle spectrum. We have therefore shown bounds 
obtained with the largest and smallest possible values of e em . Additional constraints from 
the observed Planck spectrum of the CMB have been updated recently [74]. We find that 
the new CMB limits reduce the allowed region of (mg,m f ) values only slightly when the 
conservative BBN bounds are considered. The severe BBN bounds however are more severe 
than the CMB limits in any region of the (m^,mf) plane. 

We have obtained new mass bounds from our result for ehad- By comparing our new 
bounds with the outdated ones, we find that the previous estimate of ehad from Ref. [21] 
leads to overly restrictive mass bounds. In particular, superWIMP gravitino dark matter 
scenarios, in which the dark matter density is dominated by gravitinos from charged slepton 
NLSP decays [17, 18, 21], were believed to be excluded by the severe BBN bounds [22]. 
With our update of the constraints from hadronic energy release, we find that such scenarios 
will still be viable at (rnQ,nif) = (200 GeV, 2.8 TeV) for mf L <C or (mg,m?) = 

(120 GeV, 2.2 TeV) for mg^/^, even if the severe BBN bounds are valid. With 

the conservative BBN bounds, superWIMP scenarios are also viable around (m^,nif) = 
(450 GeV, 1.2 TeV) for m fl <C m% uiil or (m^m f ) = (280 GeV, 1 TeV) for m fl ps m 5liAl . 

From the cosmological constraints, we find an upper limit on the mass of the gravitino 
LSP of (100 — 700) GeV in scenarios with a charged slepton NLSP. The wide range is 
due to the aforementioned uncertainties. We have outlined potential insights from future 
experiments that can help to reduce these uncertainties. 
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In addition to scenarios in which the NLSP sleptons freeze out with a thermal abun- 
dance, we have also considered scenarios in which a given fraction / of ^cdm comprises 
of gravitinos from NLSP decays. The electromagnetic and hadronic nucleosynthesis con- 
straints on mg and rxif have been given for / = 1, 0.1, and 0.01. The constraint from ^cdm 
is respected in such scenarios by definition. Since we consider the thermal NLSP abun- 
dance as the natural one, we have indicated explicitly the range in which the NLSP yield 
Yf required for Q^ TP = /^cdm a g rees with the estimates y^ nT i <mc i'' l i anc [ y_ mr i~ mc i.Mi ^ 

In our study of the relic gravitino density from thermal production, 0^ p , we have used 
the gauge-invariant result of Ref. [16]. While the final result for fi^ p was only given in 
the limit mg <C m g in Ref. [16], we have presented results obtained from the more general 
expression valid also for mg « m g [25]. We have shown explicitly that the contribution 
from the spin-3/2 components of the gravitino becomes important for mg > 0.5 m^. 

The dependence of on the reheating temperature after inflation Tr has allowed 
us to provide upper limits on Tr for scenarios in which a given fraction (1 — /) of ^cdm 
comprises of thermally produced gravitinos. We have described how these limits can be 
tightened by constraints from late slepton NLSP decays. Implications of the Tr limits have 
been discussed for successful thermal leptogenesis which requires Tr > 3 x 10 9 GeV [48]. 
While thermal leptogenesis cannot explain the baryon asymmetry in super WIMP gravitino 
dark matter scenarios with ^qdm ~ ^ TP > we ^ n< ^ * na * ^ ne rec i u i re d reheating tempera- 
tures become allowed for mg > 50 GeV and m g < 1 TeV in scenarios with ngg M » n? p . 

We have computed the present free-streaming velocity of gravitinos from thermal pro- 
duction and of gravitinos from NLSP decays. We find that thermally produced graviti- 
nos with mg > 100 keV can be classified as cold dark matter since they have a negli- 
gible free-streaming velocity today. In contrast, gravitinos from late decays of a charged 
slepton NLSP can be classified as warm dark matter in large regions of the parameter 
space [13, 58, 59, 27]. We have listed upper limits on the present free-streaming velocities 
from simulations and observations of cosmic structures. For superWIMP gravitino dark 
matter scenarios with ^qdm ~ ^g TP > we find that these upper limits exclude slepton NL- 
SPs with a mass of raj < 0.5 TeV, and even raj values up to about 1 TeV could be excluded 
with a better understanding of cosmological reionization and the Next Generation Space 
Telescope. However, natural superWIMP gravitino dark matter scenarios require anyhow 
mj > 1 TeV. As warm dark matter, the superWIMP gravitinos could explain the matter 
power spectrum on small scales, which is significantly lower than expected from N-body 
simulations of cold dark matter [58, 59]. 

For gravitino LSP scenarios with a charged slepton NLSP, we have outlined the possible 
interplay between astrophysical investigations and studies at future colliders. If not too 
heavy, charged NLSP sleptons will appear as long-lived charged particles in planned collider 
detectors. Such signals would point to an extremely weakly interacting LSP and could allow 
one to measure my [41]. Assuming the gravitino LSP scenario, this raj value can be used to 
refine the upper bound on mg from the cosmological constraints presented in this paper. 
If, in addition, the lifetime of the NLSP, tj, can be measured, one will be able to determine 
the mass of a possible gravitino LSP from the SUGRA prediction for tj [41, 42]. We have 
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emphasized that an experimental determination of (mg,mj) would impose upper limits on 
the gravitino density from NLSP decays, £^ TP , and thereby improve the upper limit on 
the reheating temperature after inflation Tr. 

We have described the possible identification of the gravitino as the LSP at future 
colliders. With an analysis of the NLSP decays in massive additional material around 
planned collider detectors [43, 44], distinguishing the gravitino LSP from the axino LSP 
could be possible [45, 25]. Moreover, with a kinematical determination of the LSP mass 
and the measured values for toj and rj, one can extract the Planck scale Mpi from the 
SUGRA prediction for rj. An agreement with macroscopic measurements of Mpi would be 
strong evidence for the existence of the gravitino LSP [42]. 

The various insights gained from the cosmological constraints and the prospects for 
collider phenomenology have been summarized in terms of ten viable benchmark scenarios 
with a gravitino LSP and a stau NLSP. We find three groups of scenarios: 

(i) Warm (super WIMP) gravitino dark matter from stau NLSP decays, ^qdm ~ ^ TP - 
These scenarios could resolve the small scale structure problems of cold dark matter 
but are associated with severe limits on Tr and a stau NLSP with a mass of m,f > 
1 TeV, which will be difficult to produce at the LHC. 

(ii) Mixed cold/warm gravitino dark matter with considerable contributions from ther- 
mal production and stau NLSP decays. In these scenarios, the very high reheating 
temperature of Tr > 3 x 10 9 GeV needed for successful thermal leptogenesis is possi- 
ble. The mass of the stau NLSP in these scenarios is typically m,f > 0.7 TeV so that 
it could still be discovered at the LHC. 

(iii) Cold gravitino dark matter from thermal production, f^DM ~ ^ n these scenar- 
ios, very high reheating temperatures of Tr > 10 9 GeV can still be allowed. Since 
the mass of the stau NLSP can be as small as 100 GeV, these scenarios are very 
promising for collider phenomenology and could be accessible even at the ILC. 

It remains to be seen if any of these scenarios can be verified in future experiments. In 
view of the upcoming collider experiments at the LHC, it will be worthwhile to further 
refine the cosmological constraints presented in this work. 
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